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Abstract 

We define geometric axioms for the metric and the connection of a spacetime where the gravi- 
tational influence of the connection may be interpreted as dark matter. We show how these axioms 
lead to the Einstein-Klein-Gordon equations with a cosmological constant, where the scalar field of 
the Klein-Gordon equation represents the deviation of the connection from the standard Levi-Civita 
connection on the tangent bundle and is interpreted as dark matter. 

This form of dark matter, while not quantum mechanical, gives virtually identical predictions to 
some other scalar field dark matter models, including boson stars, which others have shown to be 
compatible with the ACDM model on the cosmological scale. In addition, we quantify the already 
known fact that this scalar field dark matter, unlike the WIMP model, is automatically cold in a 
homogeneous, isotropic universe, sufficiently long after the Big Bang. 

With these motivations in mind, we show how this scalar field dark matter, which naturally forms 
dark matter density waves due to its wave nature, may cause the observed barred spiral pattern density 
waves in many disk galaxies and triaxial shapes with plausible brightness profiles in many elliptical 
galaxies. If correct, this would provide a unified explanation for spirals and bars in spiral galaxies 
and for the brightness profiles of elliptical galaxies. We compare the results of preliminary computer 
simulations with photos of actual galaxies. 

1 Overview 

There are three ideas, each of independent interest, which, if correct, would create a new connection 
between differential geometry and astronomy, very much in the tradition of general relativity's previous 
successes at describing the large-scale structure of the universe. We begin by discussing these ideas in 
general terms. 

Idea 1 : Natural geometric axioms motivate studying the Einstein-Klein-Gordon equations with a 
cosmological constant. The Klein-Gordon equation is a wave type of equation for a scalar field which 
we propose as a model for dark matter. While this geometric motivation is new, modeling dark matter 
with a scalar field satisfying the Klein-Gordon equation is not (see [ 1 8 1 for a survey of scalar field dark 
matter and boson stars). Hence, ideas 2 and 3 apply to these other works as well. 

Idea 2 : Wave types of equations for matter fields, such as the Klein-Gordon equation, naturally form 
density waves in their matter densities because of constructive and destructive interference, like waves on 
a pond, or the Maxwell equations for electromagnetic radiation. Unlike waves on a pond or the Maxwell 
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Figure 1: NGC1300 on the left, simulation on the right. The simulated image on the right results from 
running the Matlab function spiralgalaxy(l, 75000, 1, -1, 2000, 1990, 25000000, 8.7e-13, 7500, 5000, 
25000000, 10000) described in section]^ Left photo credit: Hillary Mathis/NOAO/AURA/NSF. Date: 
December 24, 2000. Telescope: Kitt Peak National Observatory's 2.1-meter telescope. Image created 
from fifteen images taken in the BVR pass-bands. 



equations, however, the group velocities of wave solutions to the Klein-Gordon equation (with positive 
"mass" term) can be anything less than the speed of light, and can be arbitrarily slow for wavelengths 
which are long enough. This allows for the possibility of gravitationally bound "blobs" of dark matter to 
form. In this paper we make some conjectures about these scalar field dark matter density waves. 

Idea 3 : Density waves in dark matter, through gravity, naturally form density waves in the regular 
baryonic matter. In the case of disk galaxies, where friction in the interstellar medium of gas and dust 
is important f9l, we exhibit examples where barred spiral density wave patterns form in the regular 
matter, as seen in figures 1, 2, 3, and 4. In the case of elliptical galaxies, where the interstellar medium 
is believed to be mostly irrelevant @|, we show how these dark matter density waves tend to produce 
triaxial ellipsoidal shapes for the regular matter with plausible brightness profiles, as seen in figures 5, 6, 
and 7. 

This paper is an attempt to put the above three ideas together in as precise a way as possible. The 
rules of the game we choose to play here are strict: define a concise set of geometric axioms, and then 
try to understand the implications of those axioms. The axioms we choose, stated in the next section, 
are more fundamental than defining an action. Instead, we declare the spacetime metric and connection 
on the tangent bundle as the fundamental objects of our universe, and then define the properties that the 
action for these two objects must have. In this paper we make the case that the simulated images in 
figures 1-6 as well as the computed brightness profile described by figure 7 are all consequences of these 
axioms. 

The resulting theory is a generalization of the vacuum Einstein equations with a cosmological con- 
stant, which already famously explains gravity, 73% of the mass of the universe as dark energy [ 14 1, the 
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Figure 2: NGC4314 on the left, simulation on the right. The simulated image on the right results from 
running theMatlab function spiralgalaxy( 1 , 75000, 1, -0.15, 2000, 1990, 25000000, 8.7e-13, 7500, 5000, 
30000000, 50000) described in section 5.4 Left photo credit: G. Fritz Benedict, Andrew Howell, Inger 
Jorgensen, David Chapell (University of Texas), Jeffery Kenney (Yale University), and Beverly J. Smith 
(CASA, University of Colorado), and NASA. Date: February 1996. Telescope: 30 inch telescope Prime 
Focus Camera, McDonald Observatory. 



accelerating expansion of the universe, black holes, and all other vacuum general relativity effects. If 
successful, this generalized theory would then, by describing dark matter, account for 95% of the mass 
of the universe [14J and explain some portion of the structures of galaxies. 

However, strictly speaking, quantum mechanics is not part of the theory we propose here, as should 
be expected since general relativity and quantum mechanics have yet to be unambiguously unified. 
Hence, our theory is clearly incomplete. This is not so bad considering that every theory known to- 
day is, in the strictest sense, incomplete. However, it is still a reasonable question to wonder if our 
theory does a good job of describing dark matter, even if it does not describe regular particulate matter. 
Hence, as a way of testing our dark matter model, we treat the remaining 4% of the mass of the universe, 
composed of particles of various kinds, in the traditional manner similar to test particles, but with mass, 
which arguably makes this theory compatible with AC DM models, as will be explained in section [3] 
The simulated images in figures 1-6 are pictures of the effect of the dark matter on the regular matter that 
we have sprinkled into the theory. Finding a less contrived way of getting regular matter into the next 
theory, while still respecting the idea of keeping our axioms as simple as possible, is an important open 
problem 

So, does the Klein-Gordon equation accurately describe dark matter and predict some observed prop- 
erties and structures of galaxies? In this paper we present evidence of this possibility by trying to un- 
derstand the effect that this model of dark matter would have on the structure of galaxies, which is a 
reasonable idea since galaxies have large components of dark matter. 

In doing so we have had to make approximations and educated guesses, so the comparisons in figures 
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Figure 3: NGC3310 on the left, simulation on the right. The simulated image on the right results from 
running the Matlab function spiralgalaxy( 1 , 75000, 1, -0.15, 2000, 1990, 100000000, 8.7e-13, 7500, 
5000, 45000000, 50000) described in section [53] Left photo credit: NASA and The Hubble Heritage 
Team (STScI/AURA). Acknowledgment: G.R. Meurer and TM. Heckman (JHU), C. Leitherer, J. Harris 
and D. Calzetti (STScI), and M. Sirianni (JHU). Dates: March 1997 and September 2000. Telescope: 
Hubble Wide Field Planetary Camera 2. 



1-6, while encouraging, should be taken in this context. Also, our "simulations" of galaxies only simulate 
the effect of the dark matter on the regular matter and hence are very primitive. Perhaps a better name 
would be "numerical experiments." However, one has to start somewhere, and it is already interesting 
that compelling patterns very much resembling actual galaxies have emerged. We will describe the 
models we have used in sections [4] [5] and [6] and the assumptions we have made. We will do our best 
to clearly label where we have had to make approximations and educated guesses, as well as rigorous 
arguments, so that readers may make their own judgments about what is presented here. 



2 Geometric Motivation 

Einstein's theory of general relativity was made possible by Gauss and Riemann who, decades before, 
began the field of mathematics now called differential geometry. Since then, advances in differential 
geometry have played a crucial role in understanding the implications of Einstein's theory. Einstein 
used differential geometry to make the qualitative statement "matter curves spacetime" precise, thereby 
showing that gravity results as a consequence of this fundamental idea. In contrast, Newton's inverse 
square law for gravity, while a great approximation in the low-field limit, has been shown to be false by 
measuring the precession of the orbit of Mercury, for example, as well as the bending of light around 
the Sun, which is twice what is predicted by Newtonian physics and exactly what is predicted by general 
relativity. Hence, understanding gravity would appear to require differential geometry. In light of this 
rich history of differential geometry playing a vital role in understanding gravity and the large scale 
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Figure 4: NGC488 on the left, simulation on the right. The simulated image on the right results from 
running the Matlab function spiralgalaxy(0. 1 , 100000, 1, -0.5, 20000, 19900, 50000000, 8.7e-15, 15000, 
5000, 82000000, 20000) described in section [5]6l Left photo credit: Johan Knapen and Nik Szymanek. 
Telescope: Jacobus Kapteyn Telescope. B, I, and H-alpha bands. 



structure of the universe, it seems reasonable, among other ideas, to look for geometric motivations for 
dark matter. 

The beginning point for our theory is to remove the assumption that the connection on the tangent 
bundle of the spacetime, an intrinsic geometric object second only to the metric in importance, is the 
standard Levi-Civita one. We compare this step to the jump from special relativity to general relativity, 
where the assumption that the metric of the spacetime is the standard flat one is removed. Our axioms 
then define the geometric properties that our action, which is now a function of the metric and the 
connection, must have. 

We note that Einstein and Cartan famously played around with ideas similar to these by removing the 
assumption that the connection was torsion free, while still assuming metric compatibility. However, as 
our beginning point, we make neither assumption. Also, Einstein and Cartan were not trying to describe 
dark matter and thus had different objectives in mind. 

Throughout this paper, the fundamental objects of our universe will be a spacetime manifold N with 
a metric g of signature ( — h ++) and a connection V. We will assume that is a smooth manifold 
which is both Hausdorff and second countable, which, while standard, deserves contemplation, as do 
all assumptions. We refer the interested reader to ll26l as an excellent reference for the fundamentals of 
differential geometry. We will also assume that g and V are smooth. These preliminary assumptions 
could be considered our "Axiom 0." 

A smooth manifold is a Hausdorff space with a complete atlas of smoothly overlapping coordinate 
charts ||26ll . Hence, we see that coordinate charts are more than convenient places to do calculations, but 
are in fact a necessary part of the definition of a smooth manifold. Given a fixed coordinate chart, let 
{d'i}, < i < 3, be the tangent vector fields to N corresponding to the standard basis vector fields of the 
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Figure 5: M87 on the left, simulation on the right. The simulated image on the right results from running 
the Matlab function ellipticalgalaxy(0. 1 , 100000, 1, 1, 2000, 1990, 50000000, 8.7e-13, 15000, 1500, 



2174670000, 1000, 2000) described in section 6.1 Note that this simulated image is the top view of 



the same simulation in the next figure which shows the side view. Left photo credit: ING Archive and 
Nik Szymanek. Date: 1995. Telescope: Jacobus Kapteyn Telescope. Instrument: JAG CCD Camera. 
Detector: Tek. Filters B, V, and R. 



coordinate chart. Let gij = g{di,dj) and Tijk = giVg^dj, dt), and let 

M = {gij} and C = {Tijk} and M' = {gij^k} and C" = {Ty-fc,/} 

be the components of the metric and the connection in the coordinate chart and all of the first derivatives 
of these components in the coordinate chart. We are now ready to state our central geometric axiom 
which motivates the remainder of this paper. 

Axiom 1 For all coordinate charts ^ : Q C N ^ R'^ and open sets U whose closure is compact and in 
the interior ofQ, [g^V) is a critical point of the functional 

F<s>,ui9, V) = /" QuadM^M' U M U C U C) dV^i (1) 
J'S'iU) 

with respect to smooth variations of the metric and connection compactly supported in U, for some fixed 
quadratic functional QuadM with coefficients in M. 

Note that we have not specified the action, only the form of the action. As is standard, we define 

Quady ({rr,}) = ^ F°^(y)x«x^ (2) 

for some functions to be a quadratic expression of the {xa} with coefficients in Y. 
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Figure 6: NGCl 132 on the left, simulation on the right. The simulated image on the right results from 
running the Matlab function ellipticalgalaxy(0.1, 100000, 1, 1, 2000, 1990, 50000000, 8.7e-13, 15000, 
1500, 2174670000, 1000, 2000) described in section [61^ Note that this simulated image is the side view 
of the same simulation in the previous figure which shows the top view. Left photo credit: NASA, ESA, 
and the Hubble Heritage (STScI/AURA)-ESA/Hubble Collaboration. Acknowledgment: M. West (ESQ, 
Chile). 



The implications of the above axiom when the connection is removed have been long understood. 
When the integrand in Axiom [T] is reduced to QuadM{M'), vacuum general relativitiy generically re- 
sults. When the integrand is reduced to QuadM{M'uM), vacuum general relativity with a cosmological 
constant generically results. Here "generically" means for a generic choice of quadratic functional, so 
that the zero quadratic function, for example, is not included in these claims. These two results were 
effectively proved by Cartan [10], Weyl B2l . and Vermeil BOl and pursued further by Lovelock ll2ni . 
The point to keep in mind is that since {g, V) must be a critical point of this functional in all coordinate 
charts, then something geometric, that is, not depending on a particular choice of coordinate chart, must 
result. Hence, if we remove the assumption that the connection is the standard Levi-Civita connection, 
the above axiom seems like a reasonable place to start. 

For organizational reasons we have placed the bulk of the geometric calculations in the appendices 
where we have provided a detailed geometric discussion of the implications of Axiom [T] In this section 
we are content to report that the Einstein-Klein-Gordon equations with a cosmological constant result 
from quadratic functional compatible with Axiom [T] More generally, we conjecture this same outcome 
for a generic choice of quadratic functional compatible with Axiom [T] Explicitly, the Einstein- Klein- 
Gordon equations with a cosmological constant (in geometrized units with the gravitational constant and 
the speed of hght set to one) are 

G + Ag = SiTfio i2—-^ l—^ + f ]g> (3) 



Y2 y 

□/ = TV (4) 
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Figure 7: The right image is the simulated elUptical galaxy image from figure 5, where the blue dots 
represent stars. If one computes the distance of each of the stars from the origin in the viewing plane, the 
histogram of those computed radii is the blue curve on the left (and is closely related to the brightness 
profile of the galaxy). The red curve on the left represents what is typically actually observed in elliptical 
galaxies according to the "i?^/^ model." This is explained in more detail in section |6] 

where G is the Einstein curvature tensor, / is the scalar field representing dark matter, A is the cosmolog- 
ical constant, and T is a new fundamental constant of nature whose value has yet to be determined. The 
other constant /xq is not a fundamental constant of nature as it can easily be absorbed into /, but simply 
is present for convenience and represents the energy density of an oscillating scalar field of magnitude 
one which is solely a function of t. It is perfectly fine to set //q equal to one, just as we have done to the 
speed of light and the gravitational constant. As a final comment, note that since 



□/ = V • V/ = -^di ( vl5l 9''d,f) = g'^ id,d,f - f ) (5) 



1 

Vlfl 

equation |4] is hyperbolic in / when the metric has signature ( — h ++)■ 

As is explained in the appendices, the scalar field represents the deviation of the connection from the 
standard Levi-Civita connection. Hence, when the scalar field is zero, the connection is the Levi-Civita 
connection. More generally, there is a formula for the connection on the tangent bundle of the spacetime 



in terms of the metric and the scalar field, which is equation 101 in the appendices. 

We must admit that we do not have a definitive idea of how the connection manifests itself physically, 
other than gravitationally, which is made explicit in the Einstein-Klein-Gordon equations. For example, 
since light rays follow null geodesies, and the geodesic equation involves the connection, it may be 
possible that light rays are affected by the connection in addition to the metric. However, one can also 
think of light rays as being along paths which are critical points of the geodesic energy functional which 
does not involve the connection. This latter view of null geodesies guarantees that these curves, once 
null, stay null. Hence, there is no guarantee that the connection affects light rays, and there is even a 
reasonable argument that it does not. The question of how the connection might be detected is a very 
interesting open problem. In this paper we will simply study the huge gravitational impact that the 
connection, through the scalar field /, has in our theory. 
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At this point an aside is warranted. We have called equation [4] the Klein-Gordon equation, which it 
is, except that the constant in the equation is not called mass as is usual, nor is Planck's constant present. 
Hence, some might want to call the above equation a modified Klein-Gordon equation to be precise. 
The point we wish to make is that there are no quantum mechanical implications being made here. We 
have "coincidentally" derived an equation which also comes up in quantum mechanics. Of course this 
is actually not much of a coincidence, since the Klein-Gordon equation is one of the simplest relativistic 
equations which one could consider 

On the other hand, this is where our discussion intersects with the fascinating works of many others 
who have studied "scalar field dark matter" and "boson stars." In the boson star case, the motivation is 
quantum mechanical, so the above scalar field / is supposed to represent the overall wave function for a 
very large number of very tiny bosons with masses on the order of 10"^'^ eV (see |[T8]| for a survey). We 
relate our constant T to the mass of the Klein-Gordon equation by noting that the Compton wavelength 
in both cases is 

A = — = — ssl3 hght years (6) 
I m 

if we take m « W~^^eV, or equivalently, T « 1/(2 light years). These other motivations are interesting 
as well but should be distinguished from the purely geometric motivations provided in this paper. For 
example, in our context it seems most natural to define / to be real valued, and we do not mean to suggest 
an interpretation of / as a probability density. At the same time, we refer the reader to the survey article 
|[T8l as well as the many works cited in that survey for many excellent discussions and ideas, most of 
which apply to this work here also, insofar as they apply to the Einstein-Klein-Gordon equations with a 
cosmological constant. Similarly, the results of this paper, from this point on, apply to boson stars and 
these other theories as well. 



3 Cosmological Predictions 

In this section we discuss the compatibility of this model for dark matter with the AC DM (Lambda 
Cold Dark Matter) model of the universe ||41ll . As our theory only hopes to explain dark matter and 
dark energy, we artificially "add in" the remaining quantum mechanical particles and electromagnetic 
radiation. Equivalently, in this section we mean to consider the standard AC DM model, but with our 
scalar field model for dark matter (which satisfies the Klein-Gordon equation) instead of the traditional 
WIMP (Weakly Interacting Massive Particles) model. 

Fortunately for us, much analysis of the cosmological features of a scalar field model of dark matter 
has already been done. Our exact case is treated in [23], referred to in that paper and elsewhere as the 
V{(j)) = (jP' case (which gives rise to the Einstein-Klein-Gordon equations). More general scalar field 
potentials are considered in that paper as well as iFTSl . ||5l, |[24l . and ||25]| . but these potentials, when they 
are even functions with positive second derivative at zero, are equivalent to the (f)^ potential in the low 
field limit. 

In |[23l (see also lIH, ll24l . ll25l ). the authors explain that scalar field dark matter has the same cos- 
mology as the the standard CDM model with WIMP dark matter particles, with the main differences only 
becoming apparent on the scale of galaxies. We refer the reader to f43] for a related discussion where the 
large scale behavior of the Einstein-Klein-Gordon equations is compared to the evolution of coUisionless 
matter It would seem reasonable that there are many more interesting problems to study in these areas, 
but the author is not expert enough to comment further. 
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Furthermore, scalar field dark matter has some advantages over the WIMP model. Scalar field dark 
matter gives a nearly flat density of dark matter in the centers of galaxies |[23l . Q, as opposed to a cusp 
of density as predicted by the WIMP model, which has not been observed. Also, numerical simulations 
show that WIMPs have a tendency to clump forming structures smaller than those observed, the smallest 
observations being dwarf galaxies [18.1 . Scalar field dark matter, on the other hand, has a natural length 
scale defined by our fundamental constant 1/T. For the group velocities of scalar field dark matter 
waves to be much less than the speed of light in order to be able to form gravitationally bound systems, 
the wavelengths of our solutions must be much larger than 1/T. 

Finally, scalar field dark matter is automatically cold, sufficiently far after the big bang, assuming 
the universe is homogeneous and isotropic. In general terms, the reason for this is the wave nature of the 
scalar field. Said another way, when two solutions to the Einstein-Klein-Gordon equation are added to 
one another, their stress-energy tensors do not add, as can be seen by the nonlinearity of the stress-energy 
tensor on the right hand side of equation |3] That is, these scalar field solutions, which can equivalently 
be thought of as wave solutions, interfere with one another both constructively and destructively. The 
result is a pressure which oscillates rapidly between being positive and negative with the average pressure 
being very close to zero, sufficiently far after the big bang, in the homogeneous, isotropic universe case. 
These statements are made precise by the following theorem, proved in Appendix [D| 

Theorem 1 Suppose that the spacetime metric is both homogeneous and isotropic, and hence is the 
Friedmann-Lemattre-Robertson-Walker metric —dt^ + a{t)'^ds% where dsf. is the constant curvature 
metric of curvature k. If f{t, x) is a real-valued solution to the Klein-Gordon equation (equation^with 
mass term T ) with a stress-energy tensor which is isotropic, then f is solely a function oft. Furthermore, 
if we let H{t) = a'{t)/a{t) be the Hubble constant (which of course is actually a function oft), and p{t) 
and P{t) be the energy density and pressure of the scalar field at each point, then 



P 



where 



and 



P 1 + e 



4T2 



(7) 



(8) 



^ — tH'{t)f{tfdt 

p{t) dt, P = / P{t) dt, H' = / , (9) 



b — a Ja b — a 

where a, b are two zeros of f (for example, two consecutive zeros). 



We comment that if we crudely approximate I-?/' (t) I (10^*^ light years) ~^ and T k, (2 light years) 
then |e| w 4 • 10^^" at the current age of the universe. We leave it to others to find better approximations 
for these values, but the point is that e is very small. Hence, this scalar field dark matter model automat- 
ically gives a cold dark matter model in the homogeneous, isotropic case, sufficiently long after the Big 
Bang. 

In contrast, in the WIMP model the particles are considered independently of one another, and so 
their stress-energy tensors do add, allowing for significant pressure terms to be present, even in the 
homogeneous, isotropic universe case. Hence, unlike the WIMP model, there is no hot version of the 
scalar field model of dark matter in the homogeneous and isotropic universe case. Since hot dark matter 
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is inconsistent with cosmological observations |[T4| . this is a nice property for a dark matter theory to 
have. 

Thus, scalar field dark matter is a serious candidate for cold dark matter. Furthermore, since it appears 
that scalar field dark matter and cold WIMP dark matter give similar predictions on the cosmological 
scale, it makes sense to go to the galactic scale to look for differences. 



We now begin our analysis which shows that scalar field dark matter satisfying the Einstein-Klein- 
Gordon equations may be related to spiral structure in disk galaxies and may predict plausible brightness 
profiles for elliptical galaxies. The cosmological constant is sufficiently small to be well approximated 
by zero on the galactic scale, so we need to study solutions to the Einstein-Klein-Gordon equations 



Studying these equations is particularly challenging because of their nonlinearities. The trivial so- 
lution to these equations is the Minkowski spacetime with zero scalar field. Also, any vacuum general 
relativity solution with zero scalar field is a solution. In addition, when the scalar field is taken to be 
complex, there exist spherically symmetric static solutions [17]. As discussed in [4J, one may think of 
these static solutions as solutions where the dispersive characteristics of the scalar field, which are mild 
for long wavelengths, are balanced with gravity. There are also real scalar field solutions analogous to 
the complex ones, but these are not quite static and are only known approximately l29l . 

In 1992, 1*341, fT6|, and then |T9l, f^Ol proposed scalar field wave equations for dark matter as a way 
of explaining the observed flatness of rotation curves lH, ll35l . |[28J for most galaxies. As described 
in the survey article ifTSl . this idea has been rediscovered many times (including by the author at the 
beginning of this work). Translated into our context, the exciting fact about the static solutions discussed 
above is that they give qualitatively good predictions for the rotation curves for galaxies. Since most 
of the matter in disk galaxies is going in very circular orbits [9], it makes sense to define the rotation 
curve of a galaxy to be the velocity of the matter of the galaxy as a function of radius. It is a striking 
fact that these curves are very flat. For example, most of the mass of the Milky Way Galaxy is going 
approximately 220 km/s in circular motion about its center (9). However, visible mass is not massive 
enough to account for these rotation curves, which is one of the motivations for the existence of dark 
matter. Hence, the prediction of fiat rotation curves is certainly intriguing, to say the least. 

However, most of the above solutions, in both the real and complex case, are actually unstable H. 
On the other hand, combinations of the above solutions (with one being the stable "ground state") yield 
dynamic spherically symmetric solutions which are stable according to numerical simulations [4} and 
give somewhat flat rotation curves. Alternatively, it is suggested in ifTSl that perhaps the dark matter 
halos, like most everything else in the universe, are rotating, and that this rotation provides additional 
stability. fTSl also points out that the bottom-up structure formation scenario of the AC DM model could 
allow, or even require, that the dark matter have angular momentum. 

The answers to the questions implied in the previous paragraph are not yet clear. What is really 
needed are careful simulations of the Einstein-Klein-Gordon equations in a perturbed cosmological set- 
ting to see what typical "blobs" of scalar field dark matter look like. These simulations could then address 



4 The Galactic Scale 
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the question of what happens when two blobs collide, whether they have enough dynamical friction to 
combine, and how they carry angular momentum. It may be necessary to simulate the regular matter 
at the same time so that energy from the dark matter can be transferred to the regular matter and then 
dissipated through friction and radiation. Also, regular matter could help stabilize the galactic potential 
which could help to stabilize the dark matter at sufficiently small radii. There are many important open 
problems in these areas. 

On the other hand, it is very easy to write down solutions to the Klein-Gordon equation if we fix the 
spacetime metric as the Minkowski spacetime. The Klein-Gordon equation then becomes 

-^ + A.)/ = T2/ (12) 

where A^, is the standard Laplacian on R^. Solutions can then be expanded in terms of spherical har- 
monics to get solutions which are linear combinations of solutions of the form 

/ = Acosicot) ■ Yn{9,(P) ■ r" • /^,„(r) (13) 

where 

/:.n(0 + ^^^/:,„(r) = (T^ - u;')f.,n (14) 
r 

and Yn{0, 0) is an nth degree spherical harmonic. Note that we require „(0) = but have not specified 
an overall normalization. Naturally, to get a complete basis of solutions we also need to include solutions 
like the one above but where cos{ujt) is replaced by sin(a;t). We will study real solutions in this paper, 
but complex solutions are quite analogous. 

It is also easy to write down solutions to the Klein-Gordon equation if we fix the spacetime metric to 
be a spherically symmetric static spacetime. For our purposes, suppose we approximate the spherically 
symmetric static spacetime metric as 

ds"^ = -V{rfdt^ + y(r)-2 {dx^ + dy"^ + dz^) (15) 

as is standard. The function V{r) acts likes the gravitational potential function from Newtonian physics. 



In this case, the above form in equation [T3]is still fine, but equation [14] is modified to become 

Vi^r? [f:,^{r) + ^^/:,„(r)) = (t^ - ^) (16) 

Since dark matter, which makes up most of the mass of most galaxies, is known to be mostly spherical 
||9l . even when the regular matter is not, the spacetime metric of a typical galaxy can be approximated 
reasonably well by one which is spherically symmetric and static. Hence, we want to understand the 
spherically symmetric static case with a potential well V{r) coming from the mass distribution of a 
galaxy as well as we can. 

There is one main qualitative difference between the Minkowski spacetime case and the spherically 
symmetric static case. In the Minkowski case, r^/^j „(r) always decays sinusoidally with amplitude 
decreasing like 1/r. Energy density decays like the square of this, which means that solutions do not 
have finite total energy. In the spherically symmetric case, there are still solutions like these, but there 
are also solutions which, if u) is small enough, have an exponential behavior as r goes to infinity. For 
most values of uj, this exponential behavior is actually increasing. However, for a discrete set of values 
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of ljj, the coefficient in front of the increasing exponential term is zero and all that is left is a decaying 
exponential behavior. The corresponding energy density is also exponentially decaying, and so solutions 
with finite total energy exist. 

The key to understanding this behavior is the term ^T^ — ■ When it is negative, the solution 

fui,n{f) is in its oscillatory domain. When this term is positive, fui,n{i") is in its exponential domain. 
The cases we are interested in are when a solution fui,n{r) begins in its oscillatory domain at r = but 
ultimately ends up in its exponential domain as r goes to infinity, with exponential decay and hence finite 
total energy. Hence, we must choose in the range 

T V{0) < < T lim V{r). (17) 

For such a solution fuj,n{'''), it is natural to define the radius Rad(/tj „(r)) to be the largest value of r for 
which ^T^ — yjyyi^ = 0. In the spherically symmetric case with positive matter density, V{r) is an 
increasing function of r, so there is a unique value for r which makes this expression zero. The point to 
take away from all of this is that for r > Ra.d{fu},n{r)), we get rapid decay for fu},n{r)- We also comment 
that for V{r) which are sufficiently asymptotic to a — 6/r at infinity, for example (the standard scenario), 
and for each n, there are an infinite number of discrete values of io which give finite energy solutions. 



with the values of these u converging to the upper limit in inequality 17 and the radii Rad(/aj,n(?')) of 
these solutions converging to infinity. We refer the reader to |[T1 for more analysis of these types of 
equations. 

Furthermore, there is an analogy between Minkowski solutions and solutions in a spherically sym- 
metric potential with finite total energy. Without loss of generality (by rescaling the spacetime coordi- 
nates), it is convenient to take V{0) = 1. Furthermore, V{r) « 1 since we are in the low field limit 



of general relativity. Hence, the range of allowable values for uj given by equation 17 is quite narrow. 
Furthermore, if we compare equations [14] and [16] we see that solutions with the same lo will start out 
very nearly the same for small r since V{0) = 1. They will begin to differ more and more as r increases, 
until finally the solution in the potential well will stop oscillating and decay rapidly to zero. From these 
qualitative observations, we make the claim that we can get a qualitative approximation for the solution 
in the spherically symmetric potential well by taking the solution in Minkowski space and modifying it 
by declaring it to be zero outside of some radius. The purpose of this approximation is not to describe the 
nuances of the solutions, but to get something that is qualitatively correct and has qualitatively similar 
properties. In this way we do not need to know the specifics of the potential of a galaxy, just the radius 
at which we wish to cutoff the matter density of the scalar field dark matter solutions in the Minkowski 
spacetime. 

At this point we keep our promise from the end of the first section and clearly state that we have made 
three approximations. First, we are approximating solutions in a spherically symmetric potential well 
with Minkowski solutions which we will arbitrarily cutoff at some radius. Second, the exact solutions 
in a spherically symmetric potential well would still be approximations in the cases we are interested in, 
since the model for rotating dark matter we are using does not give a spherically symmetric dark matter 
density or potential, although one could approximate the potential as such. Third, we have not shown 
and do not know for certain that there even exist solutions to the full Einstein-Klein-Gordon equations 
which are qualitatively similar to what we are doing here. Determining the answer to this last issue is a 
very important open problem which could benefit from simulations. 

As a next step in the near future, it will be a good idea to actually specify a spherically symmetric 



potential of a galaxy that one is trying to model, and then solve equation 16 with that potential. It will 



13 





Figure 8: Exact solution to the Klein-Gordon equation in a fixed spherically symmetric potential well 
based on the Milky Way Galaxy at t = 0, t = 10 million years, and t = 20 million years. The pictures 
show the dark matter density (in white) in the xy plane. This solution, which one can see is rotating, has 
angular momentum. 



then be desirable to match up the resulting dark matter density (which may be aspherical) with the closest 
spherically symmetric potential. This is a reasonable approach, but slightly more complicated than what 
we do here. Hence, we have made the decision to keep things as simple as possible for this preliminary 
analysis, which in this case suggests using Minkowski spacetime solutions which are arbitrarily cutoff 
at the radius of our choosing. We justify letting ourselves choose this cutoff radius since Rad(/tj,„(r)), 
defined in the spherically symmetric potential case, takes on arbitrarily large values for every value of n. 

The form of the Minkowski solution which we will study for the remainder of this paper and on 
which our simulations are based is 

f = Aq cos(a;ot)^o,o(r) + ^2 cos(w2t - 2</.) s\T?{e)r'^ f^^,2{r) (18) 

where /t^o,o(^) ^^id fuj2,2{r) satisfy equation [14] We note these solutions fall into the form of equation 



13 since both cos{2(f>) sin^(6') and sin(2i;^) sm^{9) are second degree spherical harmonics. Hence, the 
above Minkowski spacetime solution is the sum of a spherically symmetric solution (degree n = 0) and 
two degree two solutions. We have not included first degree spherical harmonics since they would not 
preserve the center of mass at the origin. We would in general expect higher degree spherical harmonics 
in galactic solutions as well, but for this paper we have chosen to focus on the case when the second 
degree spherical harmonic terms dominate. As we will see, these solutions all have 180° rotational 
symmetry with densities which rotate rigidly (which naturally is not the case when more than two distinct 
w frequencies are present). In the future, if one would want to try to model spiral galaxies with something 
other than two arms, then these other degree spherical harmonics should be studied. 



Also, initial computations suggest that the approximation in equation 18 is a reasonable qualitative 
approximation to a solution in the fixed spherically symmetric potential case. For example, the pictures 
in figure 8 show the precise scalar field dark matter densities in the xy plane at t = 0, t = 10 million 
years, and t = 20 million years in a fixed spherically symmetric potential well based on the Milky Way 
Galaxy. Like the Minkowski solution in equation [18] the solution depicted in figure 8 is also a sum of 
a spherically symmetric solution and two degree two solutions, but in a spherically symmetric potential. 
Notice how the solutions appear to have compact support, although actually the solution has simply 
decayed very rapidly outside a finite radius. Otherwise, the interference pattern is very similar to what 
we will see with the solutions in equation [18] with the most obvious difference being in the last ring or 



two, which are a bit stretched out in the radial direction compared to our model. Since all we ultimately 
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care about are the qualitative characteristics of the gravitational potential for the dark matter that we 
ultimately produce, this minor difference is something we are willing to tolerate for now. 

The reason the dark matter density is rotating in figure 8 and in our Minkowski model is most easily 
seen by making the substitution 

'0J2 - ojq" 



a 



t 



into our expression for / in equations [18] to get 

f = Ao cos{ujot)fujo,o{r) + A2 cos(a;ot - 2a) sin^(6')r^/<^2,2(r). 



(19) 



(20) 



Note that to the extent that a stays fixed in time, then the above solution does not rotate and gives a fixed 
interference pattern since both terms are oscillating in time with the same frequency loq. Hence, since 
|cj2 — wol << Wo because of inequality 17 we see that we get an dark matter interference pattern which 
is rotating according to the formula 

L02- ^0 



-t 



with period 



Tdm 



(21) 



(22) 



u;2 - ujQ 

which may be related to the pattern period of the resulting barred spiral patterns in the regular matter. 

Our next goal is to approximate the energy density fioM due to this scalar field dark matter solution 
and then to expand it in terms of spherical harmonics. From equation[20| we have that 



/ 



cos 



{coot) [Aof^,^o{r) + A2Cos{2a)sm\eyf^,,2ir)] 
+ sin(woi) [A2sin(2a)sin2(6')rV^2,2(0] ■ 



Hence, by equation 10 



Mo 



1 



-G{dt,dt) 



Svr^o 

^ (A' + |v./p)+f 



Y2 



(23) 

(24) 
(25) 
(26) 



where we have assumed a long wavelength solution in the approximation. Again, with loss of generality, 
we will always assume V{r) ^ 1, so that ~ T a;2 by inequality 17 This allows us to think of a as 
being approximately fixed in time. Hence, 



fJ-DM 



[^o/a;o,o(0 +^2cos(2a)sin2(0)r2/^2,2(r-)]' + [^2sin(2a)sin2(^)r2/c.2,2(r)]'(27) 
^0/0.0,0(0' + Al sm\eyf^,,2{r? + 2A0A2 cos(2a) sm\ey f^,^o{r)f^,^2{r) (28) 



In order to expand the above result into spherical harmonics, we need to recall that spherical harmon- 
ics defined on the unit sphere are actually restrictions of homogeneous polynomials of the same degree 
which are harmonic in B?. For example, cos(2(/i>) sin^(6') is — y^, a harmonic polynomial of degree 
2, restricted to the unit sphere and sin(2(/)) siT?{6) is 2xy restricted to the unit sphere. To summarize 
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the quick facts that we need to know, here is a quick Ust of homogeneous harmonic polynomials: degree 
zero: 1, degree one: x, y, z, degree two: — y^, 2xy, "iz^ — r^, 2x2, lyz, and then one more of degree 
four: "ihz^ — SOr^z^ + Sr'^, where = + + z^. It is easy to check that these are all harmonic in 

and hence their restrictions to the unit sphere are spherical harmonics of the same degree. 

Since cos(2a)sin2 (61) = cos(2</>-20o) sin2(6') is a rotated degree two spherical harmonic and hence 



a spherical harmonic itself (taking 0o to be fixed), the last term of equation 28 is already in the form of 
a spherical harmonic times a function of r, as we desire. The first term is as well, since it is already a 
function of r and the zeroth spherical harmonic is the constant function one. Hence, it is only the middle 
term that we need to put into the desired form. To do this, we note that 



rUin'^e = (r^-z^f = z^-2r^z' + r^ 

= A(35^4 _ 3Q^2^2 3^4) _ ifL^2(3^2 _ ^2) ^r\l), 

105^ ^ 105 ^ ^ 105 ^ ^' 

which successfully expresses this term as the sum of three terms, each of which is a function of r times 
a spherical harmonic. Hence, 

s [/o(r) + C/2(r)(3z2-r2) + [/4(r)(35/-30r2z2 + 3r^) 

+C/2(r)(r2cos(2a)sin2(0)) (29) 



fJ-o 



where 

Uo{r) = Alf^,,o{rf + ^Ayf^,,,{rf 

U2{r) = ~AyU,{r? 

U2{r) = 2ylo^2A„o(0^2,2(r-). (30) 

Our next step is to compute the gravitational potential V of this dark matter density fiDM- Notice 
that the dark matter density is the sum of a spherically symmetric term, two axially symmetric terms, and 
a term with neither of those symmetries which is rotating. As it turns out, the potential will be of this 
same form. 

We will solve for the potential by solving the Newtonian equation 

A^V = Atthdm (31) 

for each value of t, where A^^ is the Laplacian in R^. Since we are interested in solutions where the 
group velocities of our dark matter scalar field is much less than the speed of light, this Newtonian 
approximation will be fine. We will use the fact that it is easy to take the Laplacian of a function 
expressed in terms of spherical harmonics since 



(32) 
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where a is a general index and Yq. is a spherical harmonic of degree Ua- This identity follows from the 
facts that r"-"Ya{6, <j)) is a harmonic function in B?, the radial part of the Laplacian is + {2/r)dr, and 
Ax{ah) = aA^b + 2(Vj:a, V^6) + bA^a. It follows that 



V 
Anno 



Wo{r) + W2{r){3z^ - r^) + W4ir){35z^ - SOr^z^ + 3r^) 
+W'2(r)(r2 cos(2a)sin2(6l)) 



where 



with the boundary conditions 



W:{r) + ^^<(r) = Un{r) 



Wl,{0) = and lim W„(r) = 0. 



(33) 
(34) 

(35) 



Note that we mean to apply the above equation and boundary conditions to W2{r) as well. These bound- 
ary conditions produce a potential function which goes to zero at infinity. Naturally any constant may be 
added to the potential if one desires. 

Since we are going to cutoff the scalar field dark matter density outside of a radius rmax of our 
choice, for r > Vmax we will have 



<W + = 0, 

r 

which, to satisfy our boundary condition at infinity, has solution 



Wn{r) 



k 



„2n+l ' 



Hence, for r > r„ 



Wn{r) 



-<(r). 



(36) 



(37) 



(38) 



2n+ 1 

In our computer simulations, we will not actually take r to infinity, but instead, after starting with 
W^(0) = and solving the o.d.e. out to r = r^ax, add whatever constant is necessary to Wn(r) to 



satisfy equation 38 outside the dark matter radius 



Thus, in summary, given any Aq, A2,ujo,uj2 and T (which is a fixed fundamental constant of the 
universe in our theory but is not precisely known), as well as any rmax (the radius at which we cutoff 
the dark matter wave function and density), we have written down formulas in this section for everything 



else. The wave function / is defined by equations 14 and 18 The dark matter density rotates at constant 



angular speed according to equation 2 1 with period given in equation 22 The dark matter density, defined 



by equations [29] and |30j has gravitational potential defined by equations [33] and [34] Furthermore, all of 
these quantities can be computed since it is particularly straightforward to solve ordinary differential 



equations as in equations 14 andMon a computer. 



5 Spiral Galaxies 

The standard text book on galactic dynamics is the book by the same name by Binney and Tremaine 
|[9l . A great companion to this book is "Galactic Astronomy" by Binney and Merrifield lH. We also 
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refer the reader to Toomre's review article entitled "Theories of Spiral Structure" ||36]| which ||9| says "is 
still worth careful reading, even after several decades." This review article begins with three interesting 
quotes, which we repeat here: 

"Much as the discovery of these strange forms may be calculated to excite our curiosity, 
and to awaken an intense desire to learn something of the laws which give order to these 
wonderful systems, as yet, I think, we have no fair ground even for plausible conjecture." 

LordRosse (1850) 

"A beginning has been made by Jeans and other mathematicians on the dynamical prob- 
lems involved in the structure of the spirals." 

Curtis (1919) 

"Incidentally, if you are looking for a good problem ... " 

Feynman (1963) 

To these quotes we add the opening statements of Toomre's review article as a description of where 
things stood in 1977: 

"The old puzzle of the spiral arms of galaxies continues to taunt theorists. The more 
they manage to unravel it, the more obstinate seems the remaining dynamics. Right now, 
this sense of frustration seems greatest in just that part of the subject which advanced most 
impressively during the past decade - the idea of Lindblad and Lin that the grand bisymmet- 
ric spiral patterns, as in M51 and M81, are basically compression waves felt most intensely 
by the gas in the disks of those galaxies. Recent observations leave little doubt that such 
spiral "density waves" exist and indeed are fairly common, but no one still seems to know 
why. 

To confound matters, not even the A^-body experiments conducted on several large com- 
puters since the late 1960s have yet yielded any decently long-lived regular spirals. " 

Toomre (1977) 

To this 160 year old "puzzle of the spiral arms" we offer evidence of the possibility that scalar field 
dark matter density waves, which arise naturally due to the wave nature of the Klein-Gordon equation, are 
the primary driver of barred spiral density waves in disk galaxies, at least in many cases. Our proposed 
contribution to understanding disk galaxies ends there - we do not simulate the gravitational interactions 
of the stars, much less the complicated dynamics of the interstellar medium composed of gas and dust, 
nor star formation and supernova. The contribution we attempt to make here is simply to qualitatively 
model the effect of scalar field dark matter satisfying the Klein-Gordon equation on the regular matter of 
a galaxy, which we demonstrate may possibly cause barred spiral patterns in a wide range of shapes. 

If indeed dark matter is responsible for the shapes of many galaxies, this would explain why it has 
taken so long to understand spiral structure. This is especially the case if the dark matter has a dynamics 
which is different than most commonly thought. Since the bulk of dark matter has only been detected 
so far by its gravitational effects, it is reasonable to look to the details of the gravitational predictions 
made by each dark matter theory. We consider the possibility that scalar field dark matter satisfying the 
Klein-Gordon equation may have a characteristic signature which provides the beginnings of a unified 
theory for spiral and bar patterns in disk galaxies. 
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Also, we comment that we do not mean to dismiss other possible explanations of spiral structure but 
instead are simply adding another possibility which seems to hold promise. We recommend the book 
"Spiral Structure in Galaxies: A Density Wave Theory" by G. Bertin and C.C. Lin ||6l, published in 1996, 
as a fascinating description of the progress of the Lin-Shu density wave theory of spiral structure. Our 
dark matter density wave theory of spiral structure has much in common with this famous theory as both 
theories try to explain observed density waves in the visible matter. In the Lin-Shu theory, these waves 
are induced by the regular matter itself which plays the dominant role. In our dark matter density wave 
theory, we suggest that the dark matter provides the primary effect leading to barred spiral structure, for 
two reasons: there is more dark matter than regular matter, at least for large radii, and scalar field dark 
matter may have a greater tendency to produce density waves because of its wave nature. Of course in 
our theory regular matter still has gravity too, so this effect must still be considered. Hence, we see an 
opportunity for a combined theory of spiral structure coming out of our proposal which accounts for the 
gravitational effects of both the dark matter and the regular matter. Of course for a complete picture, it 
seems very likely that the interstellar medium, star formation, and supernova will need to be modeled 
too. Until all of these elements are put together into a convincing model, it will not be clear what the 
ultimate implications of the preliminary simulations we do here will be. 

We begin with a few qualitative conclusions. By looking at figure 8, for example, we see that the 
dark matter density deviates from spherical symmetry and rotates rigidly for these very special solutions. 
Thus, the level sets of the resulting potential function will not be spheres but instead will resemble 
triaxial ellipsoids. Hence, the potential function for our dark matter can be thought of as a rotating 
triaxial potential. 

Another way to get a rotating triaxial galactic potential is to have another galaxy pass nearby. In |[37l . 
Toomre shows how quadrupole forces (modeled by positive masses at and 180 degrees and negative 
masses at 90 and 270 degrees rotating around the origin), which also produce something resembling 
a rotating triaxial potential, can cause surprisingly strong spiral patterns in test particles which were 
otherwise rotating in circular motion around the origin. Hence, before we even begin our simulations, 
Toomre's work suggests that we are likely to get spiral patterns to emerge. 

As described in [9], one characteristic of disk galaxies is that they have significant amounts of gas 
and dust. Unlike stars which effectively never collide and hence do not have friction, gas and dust do 
have friction. In fact, the presumption is that it is precisely this friction which causes the collapse of 
gas and dust clouds into disks in the first place. Friction decreases total energy but conserves angular 
momentum, and a disk configuration with everything rotating circularly in the same direction allows for 
relatively large total angular momentum for its total energy since the angular momentum vectors of the 
individual masses are aligned. And sure enough, most of the visible mass of disk galaxies are indeed 
moving in roughly circular orbits. In addition, the surface brightness of disk galaxies is approximately 
modeled by 

I(i^) = V-«/^^ (39) 

where Rd is called the disk scale length of the galaxy. Hence, our simulations begin with a large number 
of point masses representing the mass of the galaxy in circular motion with area density modeled by the 
above equation. Of course since we are not modeling the gravitational effects of the regular matter, its 
initial density profile is not critically important. In fact, each of these point masses should be thought of 
as test particles which, for now, act as if they have zero mass in our simulations. 
The command line for our spiral galaxy simulations, all done in Matlab, is 

spiralgalaxy(dr, rmax^^o, Ao, >^2,Tdm, fJ-o, Rd, n-parUcies,Ttotai,dt). (40) 
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This Matlab .m file may be downloaded at http://www.math.duke.edu/faculty/bray/darkmatter/darkmatter. htmll 
which also contains the author's Matlab .m file, ellipticalgalaxy.m, for doing elliptical galaxy simulations. 
The input variables r^ax, ^o> ^2> 7dm> A'o. and Rd have already been defined. The input variable dr is 



the step size with which the o.d.e.s in equations [T4| and [34| are approximately solved. The user is allowed 
to enter either a positive or negative value for jiQ as only the absolute value is used for the actual value 
of /Uo> where a positive value denotes counterclockwise rotation for the regular matter and a negative 
value denotes clockwise rotation for the regular matter. Similarly, a positive value for Tdm gives coun- 
terclockwise rotation for the dark matter density whereas a negative value gives clockwise rotation. The 
input variable Upartides is the number of test particles used in the simulation, Ttotai is the total time of 
the simulation, and dt is the step size in time used to compute the paths of the test particles. 



As is clear from equation 14 let us define 



^ (41) 



for A; = 0, 2 to be the spatial wavelengths of our two terms. Then using the above two equations and 
equation |22j we can solve for ujq, L02, and T in terms of Aq, A2, and Tdm to get 

TT / 1 1 \ ^ 27r 
^^0 = i^{^-tj]Tdm-7^ (42) 



2 V'^l ^0/ Tdm 
TT / 1 1 \ „ 27r 



which is routine to derive. In this way we get to choose the dark matter pattern period, presumably 
roughly equal to the regular- matter pattern period, and the two wavelengths which manifests themselves 
somewhat in the rotation curve data. In theory, one could try to find best matches to known galaxies 
using our simulation and then use this to estimate T, although we have not made a careful attempt at this 
yet. 

5.1 Spiral Galaxy Simulation # 1 

As a first example, consider 

spiralgalaxy(l, 75000, 1, -1, 2000, 1990, 25000000, 8.7e - 13, 7500, 5000, 50000000, 10000) (45) 

which produces the simulated image in figure 1 as the t = 25 million years image of figure 16, rotated 
90 degrees counterclockwise. In figure 1, we compare this image to NGC 1300, a barred spiral galaxy of 
type SBbc. Notice that in this example we have |yl2| = |^o|- The Matlab code normalizes fuiQfl{r) and 
r^/a;2,2('") to have the same magnitudes for large r, so this choice maximizes both the constructive and 



destructive interferences of the two terms in equation 1 8 and hence roughly maximizes the ellipticity of 
the level sets of the potential function. We have observed in other simulations that this has the effect of 
making the bar in the simulation more like a bar and less like an oval (which we will demonstrate in our 
next example). 

Figure 9 shows the graphs of fujofl{r) and r"^ fui2,2{f) which, as one can see, have been normalized to 
have roughly the same magnitudes for large r. Naturally r'^fuj2,2{f) is the one which equals zero at the 
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Figure 9: Graphs of fugflir) and r'^fu}2,2{r) for r up to 75, 000 light years (left) and 22, 500 light years 
(right) in Spiral Galaxy Simulation # 1. Notice how the functions begin out of phase but are gradually 
becoming more in phase since they have slightly different wavelengths. 

origin. The input parameters Aq and A2 are the wavelengths of these functions in the Umit as r goes to 
infinity. 

Figure 10 shows planar cross sections of the dark matter density produced by our choice of input 
parameters in the xy, xz, and yz planes at t = 0. These densities rotate rigidly in time with period equal 
to Tdm, which is 25 million years in this example. The densities in figure 10, which decrease roughly 
like have been multiplied by in the plots to be visible for large radii. Note how the dark matter 
densities have been cutoff at r = 75, 000 light years. We point out that the discontinuity at that radius is 
not as extreme as it appears because of the factor. Also, there is dark matter density at the origin, but 
the factor suppresses this fact in figure 10. However, this density is smooth and bounded at the origin, 
as can be seen from direct calculation. These images may be compared to figure 8 which is very similar, 
except that figure 8 shows the densities without a factor of r^. 

One way to think of the densities in figure 10 is as two "blobs" of scalar field dark matter rotating 
around each other like a binary star. If we conjecture that the dark matter, in some situations, will choose 
a configuration which maximizes angular momentum (perhaps because of many collisions and mergers 
with other blobs of dark matter), then the configuration shown here suggests what these configurations 
could look like, at least qualitatively. 

Figure 1 1 shows the plots of the radial functions Wo{r), —r'^W2{r), 'ir^W4,{r), and r'^W2{r) which 
define the potential function V given by equation [33] The factors in front of the W functions are those 
relevant for z = 0, as can be seen in equation [33] so that the relative contributions of the terms can 
be judged. In this case we see that the spherically symmetric term is most dominant, followed by the 
rotating term defined by W2{r). 

Figure 12 shows the end result of all of these computations, the potential function V , graphed in 
planar cross sections. At first glance, this result is amazingly boring. After all, these images look very 
much like a generic perturbation of a spherically symmetric potential with a second degree spherical 
harmonic term, which of course they basically are. In the next two spiral galaxy simulation examples, it 
is even harder to see the perturbation, as the example shown here roughly maximizes the ellipticities of 
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the level sets of the potential. 

We are led to believe (although we have not studied this carefully) that generic perturbations of the 
potential of this type (defined appropriately) may lead to spiral patterns emerging. We refer the reader 
to [37 1 for more related discussion. Of course, it is only reasonable to consider potentials which come 
from physically plausible and common scenarios. Hence, a reasonable thought is that the main advantage 
that scalar field dark matter offers is a physically plausible and common way of achieving these triaxial 
potentials which cause spiral patterns. 

Figure 13 shows three approximate rotation curves for this simulation. We have not graphed the 
actual velocities of the test particles, although this is a good idea for a later version of this simulation. 
Instead, we have graphed v{r) = \/r\VV\ which, in the spherically symmetric case, gives the velocity 
of test particles in circular motion as a function of r. The three graphs show the value of this function 
along the x axis, yaxis, and the line y = x. We comment that we use the y = x velocity function as 
the initial velocity (in the xy plane, with zero radial component) of our test particles which start out in 
roughly circular motion. 

The interesting feature of the approximate rotation curves, of course, is how much they resemble 
the rotation curves of many galaxies, especially the way they are somewhat flat. As far as the author is 
aware, rotation curve data for NGC 1300, the galaxy in figure 1 to which we compare the results of this 
simulation, is not available. This is a common problem because the galaxies with the best pictures are 
usually ones which face us directly, whereas the ones for which the best rotation curve data can be found 
are ones at an angle to us, which allows for redshift data to be used to determine the velocities of the stars 
and gas and dust in the galaxy. Hence, we are left to speculate that, based on other galaxies for which 
there is rotation curve data ||2l, ||3l, ||35l . Il28l . that the estimates of the rotation curve in figure 13 have 
a very plausible general form. Of course our simulated rotation curves only account for the dark matter, 
so this must be taken into account. 

Another general comment worth making is that we have not tried to match the exact radius of our 
simulation to the radius of NGC 1300. The reason is that one can always scale everything in our simu- 
lation to any scale that one desires. Of course, the fundamental constant T would have to scale as well, 
but since we do not know what it is yet anyway, this is okay for now. Hence, in these first simulations 
we are simply trying to establish that there is the potential for making good fits of simulations to actual 
galaxies. 

In figures 14 - 18, we show the results of our simulations. In this example, 5000 test particles begin 
in circular motion with an initial exponentially decaying density as described by equation [39] which 
can be seen in the top left picture of figure 14. If the potential function where spherically symmetric 
(corresponding to A2 = 0), these test particles would continue in perfectly circular motion which would 
preserve the original density and image. However, since our potential is not spherically symmetric, the 
particles deviate from perfectly circular paths and respond to the aspherical gravitational influence of the 
scalar field dark matter's gravitational potential. This results in density waves in the test particles, as can 
be seen in the figures. Each frame represents 1 million years of elapse time. 

The image visualization process that we have used to generate the images in figures 14-18 is very 
important. When plotting many data points in the standard way, it is hard to see what the true density 
of data points is in regions where the data points overlap each other a lot. Instead, we have used the 
following method: Think of each data point as a disk of a small radius with uniform brightness. Define 
the overall brightness function to be the sum of the brightnesses of all of the test disks, and then graph 
this overall brightness function. In this way, the overall brightness of each frame in the figures should 
be the same. However, we have normalized each frame so that the maximum value of the brightness in 
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Figure 10: Plot of densities (in white) of dark matter times for Spiral Galaxy Simulation #1. The left 
column has a radius of 75, 000 light years and the right column is zoomed in at a radius of 22, 500 light 
years. The top row is the density in the xy plane, the middle row is the density in the xz plane, and the 
bottom row is the density in the yz plane. The densities, which roughly decrease like have been 
multiplied by to be more easily visible. Hence, the cutoff of the dark matter density at r = 75, 000 
light years is not as severe as it appears. 
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Figure 1 1 : Plots of the radial functions composing the potential function for Spiral Galaxy Simulation 
#1. The left column has a radius of 75, 000 light years and the right column is zoomed in at a radius of 
22, 500 light years. The first row is Wo{r), the second row is —r'^W2{r), the third row is 3r^VF4(r), and 
the fourth row is W2 (r) . 
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Figure 12: Graphs of the potential function in the xy plane (top row), the xz plane (middle row), and the 
yz plane (bottom row) for Spiral Galaxy Simulation #1. The first column graphs have a radius of 75, 000 
light years whereas the second column graphs have a radius of 22, 500 light years. Notice how the top 
order appearance is like V{r) = log{r), but rounded off at the origin and perturbed everywhere else to 
give triaxial ellipsoidal level sets instead of spherical ones. 
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Figure 13: Approximate rotation curves for Spiral Galaxy Simulation #1 out to a radius of 75, 000 light 
years (left) and 22, 500 light years (right). We have approximated the rotation curves with graphs of 
Y^r|Vy| (which is exactly con^ect in the spherically symmetric case) along the x axis (in blue), along 
the y axis (in red), and along y = x {in green). 

each frame corresponds to 100% white and zero brightness corresponds to black, so the overall apparent 
brightness of each frame varies slightly. 

5.2 Star Formation in Spiral Galaxies 

Probably the most striking feature of the images in figures 14 - 18 is the development of the barred 
spiral pattern which, after roughly 25 million years, somewhat resembles the appearance of NGC 1300. 
Even more generally, it is interesting that somewhat stable patterns develop at all. At a mini-conference 
at the Fetters Research Institute in Dangriga, Belize, Arlie Fetters recognized the existence of folds in 
these images, and pointed out that folds naturally have a short time stability because of their topological 
nature. Furthermore, the density goes to infinity along a fold. To be explicit, consider the mapping 
$i : i?^ — )• represented by mapping the position of the particles at time zero to their positions at time 
t. Naturally the density a.t^t{x) at time t will increase by a factor of 1/ 1 (x) | , which goes to infinity 
along a fold in the mapping. Hence, when a generic fold in our simulation develops, it will cause a curve 
of infinite density for as long as the generic fold remains. Of course one could expand this analysis to all 
of phase space for other more general situations, but there is no need to do this here since the velocity of 
our particles is a function of position at t = 0. 

Hence, this "fold argument" predicts density waves with short time stability. Moreover, along the 
fold curve the density goes to infinity. Naturally, the theoretical prediction of regions of infinite density 
could be very important for understanding the large scale process of star formation, especially since 
there is dramatically increased star formation in the arms of spiral galaxies ||9l. We note that this same 
fold argument can be generalized to apply to three dimensional space as well. Hence, we put forth the 
possibility that dark matter enforcing a "fold dynamics" upon gas and dust clouds could be a major driver, 
possibly even the primary driver, of star formation. We leave this as a very interesting question to study. 
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Figure 14: f = to f = 11 milUon y^^s for Spiral Galaxy Simulation #1. 




Figure 15: t = 12 million years to t = 23 r^^lion years for Spiral Galaxy Simulation #1. 




Figure 16: t = 24 million years to t = 35 r^^Iion years for Spiral Galaxy Simulation #1. 




Figure 17: t = 36 million years to t 



= 47 i^j^lion years for Spiral Galaxy Simulation #1. 
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Figure 18: t = 48 million years to t = 50 million years for Spiral Galaxy Simulation #1. 
5.3 Long Time Behavior 

Before moving to the next simulation, we point out that it is not entirely clear what all of the implications 
of the first simulation are. For example, is NGC 1300 in a steady-state configuration with a precise 
pattern speed, or will it evolve at some rate similar to the pictures in figures 14-18? At this point we can 
not answer this question because the simulation just described does not account for anything other than 
the gravity of the dark matter. Hence, we really can not say anything definitive on this question. 

One possibility, which we discuss here for the sake of argument, is that spiral galaxies result from at 
least two important effects, the first being the gravity of dark matter, and the second, for example, being 
friction of gas and dust and supemovae which replenish the supply of gas and dust. This second effect, 
whatever it is exactly, would be responsible for keeping most of the matter of the galaxy in the disk with 
roughly circular motion (and hence would probably not be time reversible). Ideally, this second effect 
would partially justify our choice of initial conditions in our simulations in some general sense. In this 
manner, if such a second effect like this exists, we could view a galaxy like NGC 1300 as being the result 
of two effects fighting against each other, one pushing the state toward something similar to the t = 
initial conditions, and the other pushing distributions of regular matter into spiral patterns. This could 
conceivably lead to something close to a steady state equilibrium of a trailing spiral pattern rotating at 
a constant pattern speed. If this second effect is not time reversible (like friction), then we are immune 
from the anti-spiral theorem [22|. Hence, the existence of steady state trailing spiral patterns would 
not necessarily imply the existence of steady state leading spiral patterns, which is good since they are 
typically not observed [9J. Clearly very careful simulations are needed to shed light on these questions, 
about wliich at this point we can only speculate. 

The other issue to consider is the behavior of old stars. Old stars also are mostly in the disk of the 
galaxy going in roughly circular motion, but with a higher variance than the gas and dust and younger 
stars Does dynamical friction of these older stars with the gas, dust, and other stars or some other 
process "cool" the velocities of these older stars and keep them in the galactic disk with roughly circular 
orbits? If so, then the simulations presented here might be relevant for producing spiral patterns in the 
older stars as well. An answer to the above question is needed, at a minimum, though, before one can 
make definitive conclusions. 



31 



5.4 Spiral Galaxy Simulation #2 

As a second example, consider 

spiralgalaxy(l, 75000, 1, -0.15, 2000, 1990, 25000000, 8.7e - 13, 7500, 5000, 50000000, 50000) (46) 

which produces the simulated image in figure 2 as the t = 30 million years image of figure 23. In figure 
2, we compare this image to NGC 43 14, a barred spiral galaxy of type SBa. The main difference between 
this example and the previous simulation is that A2, which controls the size of the interference pattern 
which rotates, is now only 0.15 instead of 1. Hence, the resulting potential is much more spherically 
symmetric. (The dt step size is also 50, 000 now up from 10, 000, but this makes little difference.) 
Interestingly, we notice that this smaller value for A2 results in a more oval bar than in the previous 
example. Hence, the shape of the bar, whether narrow like a bar as in the previous example or thick Uke 
an oval as in this example, is a characteristic which can be modeled by our simulations. 

The effect of lowering A2 to 0.15 can be seen in figures 19 - 22. In figure 19, notice how the 
interference pattern, which still rotates around once every 25 million years as in the previous example, 
is much more subtle now. Also, in figure 20, the spherically symmetric term represented by Wo{r) 
dominates much more than before, which is seen by looking at the values on the y axes of these graphs. 
In figure 21, the triaxiaUty of the potential function is so mild now that it is hard to even be sure that 
the level sets are not spheres from first inspection. Finally, in figure 22, the three approximate rotation 
curves, defined as before along the x axis, the y axis, and the fine y = x, are very nearly equal, another 
consequence of being close to spherical symmetry. 

Figure 23 shows the results of the test particle simulations with the computed dark matter potential 
(figure 21) rotating rigidly with a period of 25 million years. As before, the test particles begin in 
roughly circular motion but, over time, form the patterns shown. As before, there appear to be folds 
which produce relatively dense spiral arms as well as a general concentration in the central oval bar 
region. 

We use this example to make three more points: First, in our spiral galaxy simulations (unlike the 
elliptical galaxy simulations to come), we assume an exponential initial density in the test particles. We 
do not derive this density, nor do we explain it with our model. We take this initial density as given, 
presumably from other important effects like friction present in spiral galaxies which we do not model. 
Hence, the fact that the brightness goes to zero for large radii is not something that we have attempted to 
model for spiral galaxies (although we will attempt this for elliptical galaxies). 

Second, our initial distribution is not smooth but instead is composed of 5, 000 tiny bright disks. 
While this visuaUzation method works well for us, this technique also produces artificial patterns in the 
pictures which are not physical but instead reflect the initial regular pattern in which these tiny bright 
disks were originally placed. 

Third, what is real, we beheve, is that the background material not in the spiral pattern does appear 
to be a larger percentage of the test particles than in the previous example. Hence, this also appears to be 
a feature which is controlled by this model. 
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Figure 19: Spiral Galaxy Simulation # 2: Graphs of /ajo,o(^) and r^/^j 2('') for r up to 22,500 light 
years (top left). The other three images, each with a radius of 22, 500 light years, are plots of the dark 
matter density (in white) times in the xy plane (top right), in the xz plane (bottom left), and in the yz 
plane (bottom right). The densities, which roughly decrease like have been multiplied by to be 
more easily visible. Notice how the interference pattern is much subtler since | ^2/^0 1 = 0.15 now as 
opposed to 1 in the previous simulation. 
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Figure 20: Plots of the radial functions Wo{r) (top left), —r'^W2{r) (top right), 3r^W4(r) (bottom left), 
and r'^W2{r) (bottom right) composing the potential function for Spiral Galaxy Simulation #2. Notice 
how the spherically symmetric contribution given by Wo{r) dominates the rotating component described 
by r'^W2{r) which in turn dominates the remaining two terms. 



34 




Figure 21: Graphs of the potential function in the xy plane (top row), the xz plane (middle row), and the 
yz plane (bottom row) out to a radius of 22, 500 light years for Spiral Galaxy Simulation #2. The second 
column is the same as the first column except that the point of view is looking straight down so that we 
can see the level sets of the potential function in each plane. There is some distortion in the image by the 
Matlab graphics algorithm in that the domains on the right are actually perfect squares, not rectangles. 
From this we can deduce that the level sets are slightly ellipsoidal. 
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Figure 22: Approximate rotation curves for Spiral Galaxy Simulation #2 out to a radius of 22, 500 light 
years. We have approximated the rotation curves with graphs of y^r| Vy| (which is exactly correct in 
the spherically symmetric case) along the x axis (in blue), along the y axis (in red), and along y = x 
(in green). Notice how the three curves are much more similar than in the previous simulation since the 
potential is closer to being spherically symmetric. 
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Figure 23: t = to t = 50 million years (in steps ^ 5 million years) for Spiral Galaxy Simulation #2. 



5.5 Spiral Galaxy Simulation #3 

As a third example, consider 

spiralgalaxy(l, 75000, 1, -0.15, 2000, 1990, 100000000, 8.7e - 13, 7500, 5000, 95000000, 50000) 

(47) 

which produces the simulated image in figure 3 as the t = 45 million years image of figure 24, rotated 
90 degrees clockwise. In figure 3, we compare this image to NGC 3310, a "starburst" spiral galaxy of 
type Sbc which is forming stars at a rate much higher than most other galaxies. The dark matter density 
and the dark matter potential in this example are exactly the same as in the previous example. In fact, 
the only difference between this third simulation and the previous one is that the rate at which the dark 
matter rotates has been slowed to a period of 100 million years. This would appear to have the effect of 
unwinding the spiral arms compared to the previous example. 

We also wonder if more slowly rotating dark matter density waves could be correlated with starburst 
activity in galaxies where greater numbers of stars are being formed, as in the case of NGC 3310. This 
idea is consistent with the earlier "fold dynamics" suggestion that dark matter plays an important role 
in star formation. The slower the dark matter density waves move, the more time they have to act on 
the regular matter before the dark matter density wave passes through. This idea could conceivably be 
tested in part by seeing if there are certain morphology types of galaxies more commonly associated with 
starburst activity, although it is not yet definitively clear what the prediction, if any, might be. Of course, 
since there may be more than one cause of starburst activity and multiple conditions required (such as 
sufficient gas and dust), this is a delicate, but very interesting, question to study. 

We also point out that changing the period of the dark matter rotation changes the value of T, as seen 



in equation 44 In the first two simulations, T^^ ^ 10.1 light years. In this third simulation, T^^ « 2.53 
light years. On the other hand, we have not paid much attention to the overall scale of these galaxies. 
That is, we have noticed that these simulations match the shapes of certain galaxies, up to scale. To 
actually estimate the value of T would both require finding the range of values of our input parameters 
which give best fits for these galaxies, up to scale, and then scaling the resulting output morphologies to 
fit the actual galaxies. As we said before, this may be a good project, but we do not attempt this in this 
paper. 

Finally, one thing that we must say about these simulations in general is that every run, for wide 
ranges of values for the input variables, produces simulated images which look a lot like galaxies. That 
is, we have not "cherry picked" the images shown in this paper, other than to pick the ones which look as 
much as possible to example galaxies for which we happen to have good pictures. The images we have 
not picked look a lot like galaxies too. Of course, since galaxies come in such a wide variety of shapes, 
this does not prove anything. 

On the other hand, not every galactic morphology is realized by these simulations, nor was that 
expected. For example, we effectively built in a 180 degree symmetry into our model by only looking 
at a spherically symmetric component plus a second degree spherical harmonic term. The next logical 
step is to widen the class of solutions to the Klein-Gordon equation that are considered so that additional 
morphologies of galaxies can be attempted to be modeled. 
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Figure 24: t = to t = 55 million years (in steps ^ 5 million years) for Spiral Galaxy Simulation #3. 




Figure 25: t = 60 million years to t = 95 million years (in steps of 5 million years) for Spiral Galaxy 
Simulation #3. 
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Figure 26: Spiral Galaxy Simulation #4: The dark matter density times in the xy plane (top left), 
the potential function in the xy plane (top right), the level sets of the potential function in the xy plane 
(bottom left), and the rotation curve (bottom right), all to a radius of 45, 000 light years. 
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5.6 Spiral Galaxy Simulation #4 

As a final example, consider 

spiralgalaxy(0.1, 100000, 1, -0.5, 20000, 19900, 50000000, 8.7e - 15, 15000, 5000, 82000000, 20000) 



which produces the simulated image which we compare to NGC 488 in figure 4. We only include one 
other figure for this example, namely figure 26. As can be seen in that figure, the main difference in this 
example from the previous ones is that we have increased the wavelengths Ai and A2 by a factor of 10. 
As a result, the scale of the wavelengths is approaching the scale of the galaxy itself. We have not done 
enough test runs of our simulations to draw any careful conclusions about exactly what differences this 
implies, but our first impression, based on figure 4, is that this is one way to get tightly wound spirals. 

We note that with these choices of wavelengths Ai, A2 and dark matter period To a/ that T^^ ?a 512 
light years, which is much different from our first three examples. However, we can always tweak the 
values of Ai, A2, and T^m later. At this stage we are simply trying to understand the full range of the 
possible outcomes of these simulations. 

Finally, our earlier point about not optimizing our input parameters to match the photos of actual 
spiral galaxies in figures 1-4 is particularly relevant in this last example. There was a resemblance to 
NGC 488, so we put the simulated image next to this galaxy for the sake of argument. We only tried 
this one simulation in this range of values, so it is very likely that the input parameters can be adjusted, 
possibly by a lot, to create a better match. 

Each run of these simulations, when done on the author's PC at home, takes all night typically, so 
we have by no means exhausted all of the possible scenarios to try to model. Of course, converting this 
code to a faster programming language would speed things up greatly. Using multiple processors, as the 
author has recently begun doing at Duke University, is also a big help. 

6 Elliptical Galaxies 

In contrast to spiral galaxies, gas and dust are not thought to play a major role in the overall dynamics 
of elliptical galaxies |9|. Furthermore, to a good approximation in many cases, elliptical galaxies are 
mostly transparent [8]. Light from the stars of an elliptical galaxy is mostly not absorbed into the dust of 
the galaxy but instead mostly makes it out to be observed by telescopes. Hence, the observed brightness 
of elliptical galaxies is a good indiction of the true brightness. 

In addition, since gas and dust play a greatly reduced role, elliptical galaxies should be fairly well 
modeled by n body simulations. In retrospect, there is a fairly strong argument for studying elliptical 
galaxies first, before disk galaxies, as a way of testing a dark matter theory because elliptical galaxies are 
easier to model. While spiral galaxies seem more spectacular at first glance, there is plenty to understand 
about elliptical galaxies as well. 

For example, section 4.3 "Photometry of Elliptical Galaxies" of f8l is a treasure trove of data about 
elliptical galaxies. One of the most striking facts is that the surface brightness profiles of observed images 
of many elliptical galaxies, especially NGC 1700 for example, are well modeled by 



(48) 



I{R) = aexp{-5i? 




(49) 



for constants a, 6. It turns out to be more convenient to redefine the constants such that 



I{R) = /eexp{-7.67[(i?/i?e)i/^ - 1]}, 



(50) 
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where 7.67 has been chosen so that half of the total brightness is contained inside radius Re and /g is the 
brightness at this radius. That is, 

rRe poo 

2 I{R)27rRdR= I {R) 2itR dR = 7.22tt Rjle, (51) 
Jo Jo 

as explained in fE\. Naturally this model works better for some elliptical galaxies than others, and in 
some cases there are significant deviations [8J. 

We refer the reader to |[39l as an example of previous work trying to explain these observed brightness 
profiles of elliptical galaxies. Most other models depend greatly on the process by which the elliptical 
galaxy is formed, which translates into the initial conditions of the simulation. We are also faced with this 
question. However, we find that our case is qualitatively different. The dark matter density waves, which 
formed spiral patterns in disk galaxies, now effectively "randomize" the energies of the test particles in 
elliptical galaxies. Some of the test particles gain energy and are slung out to larger radii, while some of 
the test particles lose energy and spend their time at smaller radii. That is, unlike the traditional WIMP 
dark matter case where the dark matter is taken to create a roughly static potential, our rotating dynamic 
dark matter potential removes the conservation of energy constraint for each test particle. Physically, 
this can be thought of as a transfer of energy between the test particles and the dark matter. While this 
happens in any dark matter theory, the scale of this transfer appears to be much larger for our rotating 
scalar field dark matter model than for the WIMP model. 

To see if this model for rotating scalar field dark matter could reproduce something close to the 
observed brightness profiles, we created a Matlab simulation whose command line is 

ellipticalgalaxy((ir, rmax,Ao, A2, Xq, A2, Tda/, /Uq, Ri, Upartides^Ttotah dt, sun)- (52) 

This Matlab .m file may be downloaded at"http://www.math.duke.edu/faculty/bray/darkmatter/darkmatter.html 
which also contains the author's Matlab .m file, spiralgalaxy.m, for doing spiral galaxy simulations. 

Notice that almost all of the input variables are the same as spiralgalaxy.m. In fact, the model for the 
rotating dark matter and the corresponding rotating dark matter potential are exactly the same. The only 
differences, really, are the initial conditions that we choose for our point particles and how we display 
the results. 

Notice that we no longer have Rd as an input parameter, which was the disk scale length of the spiral 
galaxy, which determined the initial density of the test particles. Instead we are trying to determine 
equilibrium densities for test particles in mostly radial motion, as is the case for elliptical galaxies [8]. 

This raises the question of what the initial placement of the test particles should be. We chose 
something very simple - place all of the test particles on a sphere of radius Ri, the initial radius, but 
with zero velocity. When the test particles collapse down toward the origin, the result will be a kind of 
explosion as the stars come back out. Elliptical galaxies are thought to be the result, at least in many 
cases, of other galaxies which have merged, perhaps many times t9J. These initial conditions are inspired 
by this as well. 

The other input variable is sun which has nothing to do with the actual simulation but only with 
how the results of the simulation are displayed. Since we want to understand the predicted brightness 
profile, we will want to know the distribution of the test particles in the viewing plane. To do this, we 
determine the radius of each test particle in the viewing plane and then create a histogram of the resulting 
radii, where the bin size for each interval of radii is s^n- We will refer to these histograms as graphs of 
the "radial frequency in the xy plane." We will then compare our measured radial frequencies with that 
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Figure 27: The initial positions of the test particles for Elliptical Galaxy Simulation #1 are on a sphere 
of radius 15, 000 light years. The right image is simply a blurry version of the left image. The initial 
velocities are all zero. 



which is predicted by equation 50 which we will approximate by 

radial frequency I{R) • IttR ■ sun- (53) 

This is easily done by defining to be the median value of the measured radii in the viewing plane and 
by choosing le to give the correct overall normalization so that the total brightnesses are the same. 

We remind the reader that the test particles in these simulations, which are mostly supposed to repre- 
sent stars in this case, are still treated as though they have zero mass. Hence, we are only measuring the 
gravitational effect of the dark matter on the stars. This is instructive for a first simulation. It will also 
be roughly accurate in cases where the dark matter makes up most of the mass of the galaxy. However, 
simulations which do account for gravitational interactions between stars is a logical next step. 



6.1 Elliptical Galaxy Simulation #1 

We will focus on a single elliptical galaxy example in this paper which was used to create the simulated 
images in figures 5 and 6 using the following command line in Matlab: 

ellipticalgalaxy(.l, 100000, 1, 1, 2000, 1990, 50000000, 8.7e - 13, 15000, 1500, le20, 1000, 2000). 

(54) 

We refer the reader to Spiral Galaxy Simulation #1 and figures 9, 10, 11, 12 for qualitatively correct 
pictures of the dark matter density and its corresponding potential which will be almost exactly the same 
as in this example (except that we made a minor change in the dark matter radius from 75, 000 light years 
to 100, 000 light years). The sign of A2 being switched simply rotates the potential by 90 degrees in the 
xy plane, which is also irrelevant for our purposes. The only other difference is that in this case the dark 
matter potential will rotate slower with a dark matter period of 50, 000, 000 million years. 

As previously mentioned, the initial placement of the test particles is on spheres, which in this case 
will have a radius 15, 000 light years as seen in figure 27. These spheres naturally collapse toward the 
origin and then return back to something close to their original positions, but modified because of the 
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Figure 28: The positions of the test particles for Elliptical Galaxy Simulation #1 after 2.175 billion years. 
The top row is a top view and the bottom row is a side view. The right column is simply a blurry version 
of the left column. 

rotating dark matter potential in the background. Over time, as previously described, a wide spread in 
the energies of the test particles develops. After a bilUon years or so, the distribution of test particles 
settles down to a kind of limiting distribution which looks qualitatively similar to elUptical galaxies, as 
can be seen in figure 28. 

We note that the images in figure 28, while very intriguing, are not perfect. The ratio of the major 
axis to the minor axis of the top view is very close to one, but this ratio is about 4: 1 in the side view 
(according to a formula used to calculate this in the simulation which weights large radii more than small 
radii). Thus, the side image appears to be too flat, since this ratio is generally not observed to exceed 
3:1. Of course our simulation is very rough, not even accounting for star - star gravitational interactions, 
so we were only hoping to get something in the right ballpark. Also, our initial conditions could be 
improved to be more reflective of how elliptical galaxies are actually formed. And finally, our model for 
rotating scalar field dark matter is only an approximation. 

Also, the side view in figure 28 is the view which maximizes the apparent ellipticity. A random point 
of view will rarely give the side view, so a typical observed ellipticity will be much reduced. We note 
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Figure 29: Elliptical Galaxy Simulation #1: The left image is the radial frequency in the viewing plane 
(in blue) compared to the standard model (in red) from equation 53 The right image is the location of 
test particles in the xy viewing plane on which the the radial frequency is based. The median radial value 
for the test particles is approximately 28, 000 light years. The arrow in the right image is irrelevant and 
simply shows the direction of the minor axis of the potential in the xy plane, which rotates. This figure 
is a repeat of figure 7. 



that a histogram of the observed ellipticities of 2, 135 elliptical galaxies is found in figure 4.33 of fS^. It 
is a logical next step to study the predicted histogram of ellipticities of any theory, so we list this as a 
good problem to study. 

We can be even more quantitative about understanding these results if we look at the radial frequency 
in the viewing plane, which is related to the observed brightness profiles of elliptical galaxies. In figure 
29, we graph the radial frequency in the xy viewing plane correspond to the top view of the simulation 
in the top row of figure 28. We find such a close match, with such a generic initial condition, to be very 
encouraging. Hence, we are very intrigued by the possibility that scalar field dark matter with angular 
momentum could be a major factor in determining the brightness profiles of elliptical galaxies. 



6.2 More Brightness Profile Data 

We comment that we have found even better matches than the one in figure 29 for the radial frequency 
of the test particles, which of course is closely related to the brightness profile of a galaxy by equation 



53 We do not go into detail about these simulations because of space considerations and because, in 
retrospect, they use an initial condition which is suspect. Namely, the test particles begin on a very small 
sphere of radius only 5, 000 light years. As a result, the test particles apparently pick up a lot of angular 
momentum and rotate about the origin more than desired. This appears to result in eccentricities higher 
than ideal. But we decided this data was worth sharing anyway, realizing that it also is not quite perfect. 
What is interesting is that many of the test particles get energized to radii many times their initial 5, 000 
light year radius well exceeding 100, 000 light years or more, as long as the dark matter has enough 
angular momentum. 

We begin by noting that if the dark matter is not rotating at all {A2 = 0), then by conservation of 
energy in a fixed potential, the (massless) test particles will never get more than 5, 000 light years from 
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Figure 30: These six radial frequency plots result from running the Matlab function described in the 
text where A2 ranges over the values 0.5 (top left), 0.75 (top right), 1 (middle left), 1.25 (middle right), 
1.5 (bottom left), and 2 (bottom right). In that same order, the radial upper limits in these six plots are, 
respectively, 40000, 220000, 400000, 500000, 400000, and 600000 Ught years. 
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the origin. If the dark matter is only rotating a little bit (|j42/j4i| << 1), not much diffusion of energies 
occurs. 

The six radial frequency plots in figure 30 result from running the Matlab function 

ellipticalgalaxy(0.1, 100000, 1, ^2, 2000, 1990, 50000000, 8.7e-13, 5000, 1500, 1000000000, 1000, 2000) 

(55) 

where A2 ranges over the values {0.5, 0.75, 1, 1.25, 1.5, 2} from the top left to the bottom right. Note 
that Ai = 1. The greatest interference pattern occurs (as a percentage of the overall pattern) when 
1^2! = |^i|> and interestingly, the corresponding radial frequency plot in figure 30 is quite a good fit. 

The fact that we get the best fit (in this imperfect example) when the interference pattern is maxi- 
mized suggests that elliptical galaxies are best modeled in this theory by dark matter with large angular 
momentum (|yl2/^i| ~ 1)- On the other hand, spiral galaxies were successfully modeled with |yl2/^i| 
both equal to 1.0 (Spiral Galaxy Simulation #1) and 0.15 (Spiral Galaxy Simulations #2 and #3). It could 
be that practically all elliptical galaxies have dark matter with large angular momentum due to a major 
collision with another galaxy in their past, whereas spiral galaxies are more likely to take on the whole 
range of angular momenta. More study is required to address this question. 

6.3 Ripples in Elliptical Galaxies 

We could not end a section on elliptical galaxies without pointing out some elliptical galaxies actually 
have visible "ripples" in their images. In fact, about 10 to 20 percent of early-type galaxies have visible 
sharp steps in their luminosity profiles, as can be seen in NGC 3923 in figure 4.40 of El- This can not 
be fully appreciated without looking at the image, so we highly suggest the reader get a copy of 1 8 1 , look 
at the image, and read the discussion. The striking characteristic of these ripples is that they tend to be 
smooth arcs of ellipses centered around the center of the galaxy. In NGC 3923, these arcs range in radii 
from 3 to 95 kpc |[8l. The ripples are on the order of 3% to 5% variations in the surface brightness of the 
galaxy. 

The issue of these ripples is not settled. Since our dark matter density has ripples in it, as seen in 
figures 10 and 19 and especially figure 26 where the wavelengths Ai and A2 have been increased, it is 
only natural to wonder if there is a connection. We leave this as a very interesting open question. 



7 Open Problems 

As exciting a project as this has been for the author, there are vastly more questions left open than answers 
to problems solved. Hence, it seems appropriate to end this paper with a list of open problems. 

Geometric Analysis Problems 

Problem 1 Prove or disprove conjecture^in appendix^ that a generic quadratic functional compatible 
with Axiom^leads to the Einstein-Klein-Gordon equations. 



Problem 2 Prove or disprove that all actions compatible with Axiom^are of the form in equation 108 



Problem 3 Classify all actions of the form in equation 108 and describe the resulting Euler-Lagrange 
equations. What kinds of phenomena exist for the solutions to these equations? 
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Problem 4 What do generic stable solutions to the Einstein-Klein-Gordon equations look like? Are there 
solutions corresponding to a gravitationally bound "blob" of scalar field with angular momentum? 



Astrophysics Problems 



Problem 5 Does dark matter enforcing a 'fold dynamics" (discussed in section 5.2) upon gas and dust 
clouds drive star formation ? 

Problem 6 If dark matter is connected to star formation, what characteristics of the dark matter are 
related to starburst activity in galaxies? Do the magnitude and speed of the dark matter density wave 
affect the star formation rate? If the magnitude and speed of the dark matter density wave can be 
associated with certain morphologies of spiral galaxies, there could be a testable hypothesis. 

Problem 7 Is there a connection between the observed ripples in early-type elliptical galaxies ( discussed 
6.3 \ and the ripples in the scalar field dark matter density? 



m section 



Problem 8 What phenomena does this scalar field dark matter predict that can be tested or detected by 
gravitational lensing? 

Problem 9 Study and explain the relationships between the input parameters of the spiral galaxy simu- 
lation in this paper and the morphologies of the disk galaxies created. 

Problem 10 Estimate the value o/T (the constant in the Klein-Gordon equation) by making best fits of 
rotation curves, brightness profiles, observed distribution of ellipticities, and spiral patterns of galaxies. 
Can the numbers and sizes of dwarf galaxies be used to estimate T? 

Problem 11 Compare the distribution of ellipticities predicted by simulations of the scalar field dark 
matter theory (from random viewpoints) to the observed distribution of ellipticities 0/2, 135 elliptical 
galaxies in figure 4.33 of [8]. Be careful to use the same procedure for defining ellipticity in the simula- 
tions as in the observations. 



Computer Simulations 

Note: http://www.math.duke.edu/faculty/bray/darkmatter/darkmatter.html has the Matlab simulations 
described in this paper available for download. 

Problem 12 Improve the elliptical galaxy simulation discussed in this paper by adding in the gravita- 
tional infiuence of the stars on each other. 

Problem 13 Generalize the simulations in this paper by incorporating higher degree spherical harmon- 
ics into the scalar field dark matter solutions. As long as no more than two frequencies are used, the 
solutions will also approximately rotate and produce potentials that approximately rotate. See if spiral 
galaxies with more than two arms can be simulated with these new rotating potentials. 

Problem 14 What are the dynamics of the stars and the gas and dust of a galaxy with a typical stable 
solution (from problem to the Einstein-Klein-Gordon equations in the background, and more specifi- 
cally, does spiral structure result for a wide range of conditions? How about observed elliptical galaxy 
characteristics ? 
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Problem 15 Does dynamical friction manifest itself in solutions to the Einstein-Klein-Gordon equations, 
and if it does, how does it do it? Is it necessary to have gas and dust around to dissipate energy to allow 
a dynamical friction process to form larger and larger gravitationally bound "blobs" of scalar field dark 
matter? 

Problem 16 Create careful simulations of the Einstein-Klein-Gordon equations in a perturbed cosmo- 
logical setting to see what typical "blobs" of scalar field dark matter look like. These simulations could 
then address the question of what happens when two blobs collide, whether they have enough dynamical 
friction to combine, and how they carry angular momentum. It may be necessary to simulate the regular 
matter at the same time so that energy from the dark matter can be transferred to the regular matter 
and then dissipated through friction and radiation. Also, regular matter could help stabilize the galactic 
potential which could help to stabilize the dark matter at sufficiently small radii. 

Physics Problems 

Problem 17 Is there a way to modify Axiom^or Axiom^to include quantum mechanics? 

Problem 18 What are the physical interpretations of the matter fields defined by Axiom^( discussed at 
the end of appendix^? 

Problem 19 How does the connection on the tangent bundle of a spacetime manifest itself physically? 

Problem 20 What is the behavior of scalar field dark matter right after the Big Bang? Note that P /p 
(from theorem^in section^ is not necessarily close to zero in this case. 

Problem 21 What are the physical consequences of the oscillating pressure of scalar field dark matter 
Note that in the real scalar field case, the pressure necessarily oscillates between being positive and 
negative, while averaging to something very close to zero, as described by theorem^in section^ 

Problem 22 How does a scalar field model of dark matter affect the black holes at the centers of most 
galaxies? At what rate does the scalar field dark matter fall into the black hole? 
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A A Brief Introduction to General Relativity 



In this appendix we provide a brief introduction to general relativity relevant to the next section where 
we show how the Einstein-Klein-Gordon equations with a cosmological constant result from Axiom [T] 
The differential geometry discussed in these two appendices is detailed but is also very standard. Our 
discussion will not require much knowledge beyond the first three chapters of |[26l . which we recommend 
as an excellent text book for those who want to use differential geometry to study the large scale structure 
of the universe. 

Special relativity, Einstein's theory which unified space and time, is really all about the study of 
the geometry of the flat Minkowski spacetime ll26l . The defining idea leading to general relativity is 
the idea that spacetime is not necessarily flat. Philosophically, the question "Why should spacetime be 
curved?" has a great response: "Why should spacetime be flat?" What is remarkable is that removing 
the assumption that spacetime is flat leads to general relativity as arguably the simplest way of defining 
what the spacetime metric should be, as we will explain. Moreover, experiments have shown that general 
relativity describes gravity as accurately as can be measured. 

The central equation of general relativity is the Einstein equation 

G = 8ttT, 

where G = Ric — ^Rg is the Einstein curvature tensor and T, called the stress energy tensor, describes 
the local energy and momentum density. However, it is important to realize that this equation by itself is 
only half a theory in that this equation describes how the matter curves the metric but does not explain 
how the matter evolves forward in time. From this point on we will only discuss complete theories where 
matter evolution equations are specified, fully realizing that one can always add in an unspecified stress 
energy tensor later if desired. For example, 

G = 0, 

the vacuum Einstein equation, is a complete theory in our terminology because there is no matter term 
to evolve. 

It is worth pointing out that vacuum general relativity is already an intricate and very interesting 
theory. For example, black holes are solutions to the vacuum Einstein equation when all of the matter 
has already fallen into the black holes. This demonstrates an interesting characteristic of general relativity 
- zero local energy and momentum density does not imply that the system has zero total mass. In fact, by 
the positive mass theorem [30,1 , 1311 , an isolated asymptotically flat vacuum solution must have strictly 
positive total mass, unless the spacetime is flat everywhere (and hence is the Minkowski spacetime). 

Also, matter can be approximated by black holes of various sizes. For example, one could replace 
the Sun and the planets by black holes of equal mass, and the corresponding n-body problem with black 
holes is a solution to G = 0. For that matter, one could replace every star, planet, meteor, and piece 
of dust in a galaxy with a black hole of equal mass, and the resulting system would be a solution to 
G = 0. The dynamics would be a good approximation to the original system except when the objects 
got close to one another. All of the above statements are true assuming that an existence theory for the 
vacuum Einstein equation G = generically exists, which is currently an open problem. The difficulty 
in even proving an existence theory for G = is another indication of the highly nontrivial nature of the 
equation. 

Einstein chose G = Ric — ^Rg once he learned that the divergence of G equaled zero by the second 
Bianchi identity |f26l. This zero divergence property can be interpreted as a local conservation property 
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for energy and momentum density, a understandably desirable property. Hilbert realized that the zero 
divergence property of G was a direct result of the fact that metrics which are critical points of the 
functional 

Hu{g)= [ RdV, (56) 



JU 

satisfy the corresponding Euler-Lagrange equation G = on the spacetime N. The role of the open set 
U in the above equation, whose closure is required to be compact, is to keep the integral finite. We define 
a metric (7 to be a critical point of Hu{g) if this quantity does not change to first order for all smooth 
variations of the metric with compact support in the interior of U. Note that we are not minimizing or 
maximizing anything but are simply looking for critical points. 

Moreover, Hilbert realized that the zero divergence property is true for any Euler-Lagrange equation 
resulting from an action which is invariant under reparameterizations of the metric, as the Einstein- 



Hilbert action in equation 56 is. Today we think of this fact as a special case of Noether's theorem. 
Hence, the zero divergence property of G, while nice, is nothing particularly special and is not the 
fundamental point. Instead, vacuum general relativity is best understood as defining the metric to be a 



critical point of the Einstein-Hilbert action in equation 56 This raises the truly relevant question: What 
is so special about the Einstein-Hilbert action? 

We begin by noting that in some ways the Einstein-Hilbert action is the simplest nontrivial scalar 
valued functional of a metric. The scalar curvature of a spacetime metric is a scalar at each point, so it 
is perfectly natural to integrate it. However, there is another characterizing property of this functional 
which is seen by expressing it in a coordinate chart. 

Given a coordinate chart mapping $ : — )• i?^ = { {x^ ^)}, define 

gij=g{di,dj), where di = D^-^ {J-}j (57) 

to be the components of the metric in this coordinate chart. Then standard calculations show that the 
formula for the scalar curvature in terms of the metric in a coordinate chart is 

R = - 9'^9^^)9ijM+ diiMdah^c- 

^^g^-g^^g^- - ]^g''^g^'g^^ - g^'^g^^g^^ - \g'^g'^''g^' + g'^g'^^g^'^ (58) 

where commas denote coordinate chart derivatives and the Einstein summation convention of summing 
over indices which are both raised and lowered is in effect. As usual, we define [g"^^] = [gij\~^ as 
matrices. Also, the formula for the metric volume form in terms of the metric in a coordinate chart is 

dV = \g\^'^ dVRA (59) 

where \g\ = \det{[gij\)\. Using \g\^i = \g\ g'^'^gmn,i and g'^ = -g"^g^''gab,k we can integrate by parts 
over an open set U (whose closure is compact) in the coordinate chart to get 

RdV = / {g'^9^^ - 9"^9^^)9ij,kyi\9\^'^dARA+ 9^j,k9ab,c■ 
JU Jd{^{U)) J^{U) 

. Igia^ife^fcc ^ ^giagjcgkb ^ ^gijgabgkc _ ^gijgacgkb\ |^|l/2 ^y^^ (go) 
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where f/ = z^)ij4 and v is the outward unit normal to ^{U) in the coordinate chart. 
Note that our formula for the Einstein-Hilbert action is of the form 

/ RdV= / QuadA^(M') (iVR4 + boundary term, (61) 

where 

M = {gij] and M' = {5^- J 

are the components of the metric in the coordinate chart and all of the first derivatives of these compo- 
nents in the coordinate chart, and where 

Quady ({x,}) = F'"^{Y)xaXis (62) 

for some functions {F°''^} is a quadratic expression of the {xa} with coefficients in Y. What is interest- 
ing about this particular quadratic expression is that while its value at each point depends on the choice 



of coordinate chart, by equation 61 its integral over ^{U) is invariant with respect to smooth variations 
of <I> compactly supported in the interior of U (which hence do not affect the boundary term). Even 
more intriguing, this quadratic expression is the unique one with this property, up to multiplication by a 
constant, which follows from the works of Cartan ifTOl . Weyl ll42l . and Vermeil f40] (and pursued further 
by Lovelock [21]). 
Hence, if we define 

F^,u{9) = f QuadM(M') dVn^ (63) 
for that same quadratic expression, then 

Hu{g) = F'S>,u{g) + boundary term (64) 



by equation 6 1 Hence, the Euler-Lagrange equation for F$ (g) , for variations away from the boundary, 
must be the same as the one for Hu{g), which is G = 0. 
More generally, let 

F^,u{9) = I QuadA/(M') dV^, (65) 

for any nontrivial quadratic expression, but require that gf is at a critical point of F^jj{g) for all coordinate 
charts $ and all open sets U (whose closure is compact and in the interior of O). Then the resulting Euler- 
Lagrange equation can not depend on any one particular choice of coordinate chart, and so by uniqueness 
as proved by Cartan [ 10|, Weyl [,42 J . and Vermeil Ii40,l up to a multiplicative constant, we must have 

c- Hu{g) = Fq,^u{g) + boundary term (66) 

for some nonzero c. Hence, again the Euler-Lagrange equation for F^^u{g), for variations away from 
the boundary, must be the same as the one for Hu{g), which again is G = 0. 
Now suppose we let 

F^,u{9)= [ QnadMiM' \J M) dVj^4 (67) 
J<S>{U) 

for any quadratic expression (which is nontrivial on the M' — M' part), and continue to require that g is 
at a critical point of F^^u{g) for all coordinate charts <I> and all open sets U (whose closure is compact 
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and in the interior of Then the resulting Euler-Lagrange equation still can not depend on any one 
particular choice of coordinate chart. Furthermore, the uniqueness theorems of Cartan UTOj . Weyl ||42]| . 
and Vermeil HOl imply the even more general result that 

c / {R-2k)dV= F^^u{g) + boundary term (68) 
Ju 

for two constants c 7^ and A. Hence, the Euler-Lagrange equation for F^^ij{g), for variations away 
from the boundary, will be the same as the expression on the left. The expression on the left is called the 
Einstein-Hilbert action with cosmological constant and is well known to have Euler-Lagrange equation 

G + Ag = 0, (69) 

which is the equation for vacuum general relativity with a cosmological constant A. 

When the cosmological constant term is moved to the other side of the equation and interpreted as 
matter as in G = —Ag = SirT, this matter is given the name "dark energy" (not to be confused with 
dark matter) and has the same stress energy tensor as a fluid with positive energy density and negative 
pressure (for A > 0). The idea that a fluid can have negative pressure is confusing to many, but the point 
is that dark energy is not a fluid. Dark energy is simply the phenomenon resulting from the constant in 



equation 69 above. Standard theory today assigns a very small but positive value for A to explain the 
accelerating expansion of the universe. While its density is very small, it is everywhere, and so ends up 
accounting for roughly 73% of the mass of the universe in the standard AC DM model fT4l . 

The above discussion provides the motivation for Axiom [T] from section [2] It is also helpful for 
understanding the implications of Axiom[T] which we discuss in the next appendix. 



B Derivation of the Einstein-Klein-Gordon Equations with a Cosmologi- 
cal Constant from Axiom [H. 

We begin by noting that the connection V on the tangent bundle of a spacetime iV is a fundamental 
geometric object, arguably second only to the metric g in importance. While the metric is used to 
measure length, area, volume, etc., the connection is used to differentiate a vector field in the direction 
of another vector field. The curvature tensors are then usually defined in terms of the connection, for 
example. 

Given any metric, there exists a unique standard connection called the Levi-Civita connection V 
which, in addition to the linearity and Leibniz properties that all connections satisfy, is both metric 
compatible and torsion free. It is a short exercise to show that given a general connection V, its difference 
with the standard connection (or, for that matter, any other connection) 

D{X, Y, Z) = {VxY, Z) - {VxY, Z) (70) 

is a tensor, meaning that D is multilinear in each of its three slots. Note that we have used the convenient 
notation g{X, Y) = {X, Y). We refer the reader to ||26| for more about connections. 

Given a fixed coordinate chart, let {Sj}, < z < 3, be the tangent vector fields to N corresponding 
to the standai^d basis vector fields of the coordinate chart. Let gij = g{di, dj) and Tijk = (Vd^dj^dk). 
Then the above tensor equation, in coordinates, becomes 

Dijk ^ijk '^ijk- C^l) 
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Furthermore, the Koszul formula for a connection 11261 impUes that 

f jjfc = 2 (9ik,j + 9jk,i — 9ij,k) ) (72) 

which we note has first derivatives of the metric in its expression (which will be important later). 
Axiom [T] directs us to look for functionals of the form 



where 



((,, V) = / QuadA^(M' UMUC'UC) dVnA (73) 



M = {gij} and C = {Ti,k} and M' = {gi,,k} and C = {Tijk,i} (74) 

are the components of the metric and the connection in the coordinate chart and all of the first derivatives 
of these components in the coordinate chart. From the previous appendix we know that 

i^*,c/(ff,V)= / {R-2A)dV (75) 
Ju 

falls into this form, at least up to a boundary term which we do not write down since it is irrelevant for 
the ultimate Euler-Lagrange equations which will be produced. Also, 

F<s,,u{g, V) = / {R-2A- Quad JL>)) dV (76) 
Ju 



will be of the correct form as well by equations 7 1 and 72 again, up to the usual boundary term. Here 
Qua.dg{D) means some quadratic functional of D with coefficients in the metric g. Note that the above 
two expressions have no dependence on <I>. 

However, Axiom [T] does not allow us to add in a general expression of the form Quad^(VZ)) to the 



integrand. This is because the F terms of equation 71 already contain first derivatives of the metric, 
so taking another derivative and squaring would introduce quadratic terms of second derivatives of the 
metric, which is not allowed by the axioms. However, the F terms go away if we look at the torsion part 
of the D tensor defined to be the D tensor antisymmetrized in the first two indices. The torsion tensor 
(with the third index lowered) has components 

Tijk — Dijj^ DjiP; (77) 

= i^ijk — ^ijk) — i^jik — ^jik) (78) 
— ^ijk Tjifc (79) 

by the symmetry of the f in the first two indices. 

Again, we would like to consider some kind of derivative of this torsion tensor which we could then 
square. However, Levi-Civita covariant derivatives generically would introduce gT terms which, when 
squared, would not be allowed by our axioms. Similarly, taking a V covariant derivative would introduce 
F^ terms before we squared, which is also not allowed by the axioms. As before, we need to find a way 
to take a derivative of some part of the torsion tensor which does not require a covariant derivative. 

Define 

^ijk = -^{Tijk + Tjki + Tkij) (80) 
= ^iDijk — Djik + Djki — Dkji + Dkij — Dikj) (81) 

— g ij^ijk Tjifc ~l~ -Tj fci ^kji ~l~ ^kij ^ikj ) (82) 
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to be the fully antisymmetric part of the difference tensor D (and half the fully antisymmetric part of 
the torsion tensor T. Thus, 74^^ are the components of a three form. Hence, we can take the exterior 
derivative of 7 to get 

djijkl = Ijkl,! — IkliJ + lUj,k — lijk,l (83) 

which are the antisymmetric coefficients of the tensor ^7 which do not involve derivatives of the metric, 
just derivatives of T. Hence, functional of the form 

F^,u{g,y) = ^(cii - 2A - g|d7|2 -Quadg(I))) dV, (84) 

are allowed by the axioms, up to a boundary term which is irrelevant for the Euler-Lagrange equations 
produced. 

Conjecture 1 All junctionals for which there exists a smooth metric and smooth connection which satisfy 
Axiom^with that functional are of the form of equation 84 



We leave this as a very interesting geometric conjecture to study. We will take c 7^ so that the resulting 
Euler-Lagrange equations for the metric are nondegenerate. Without loss of generality, we may as well 
take c = 1. 

We comment that for generic quadratic expressions Quadg(-D), all of the components of D other 
than the fully antisymmetric part 7 (and irreducible components of the same type) will have to be zero to 
satisfy their zeroth order Euler-Lagrange equations. The irreducible components of the difference tensor 
D of the same type as 7 are the three traces of D. For our purposes here, we will consider the case where 
the quadratic expression does not have any cross terms so that these three traces are zero as well. We 
comment that two of these traces are in the metric compatibility part of the difference tensor and one is in 
the torsion part of the difference tensor. We consider this case because it is representative of the general 
case in that a generic quadratic expression will also lead to the Einstein-Klein-Gordon equation with a 
cosmological constant. 

Hence, in this representative case, 

Dijk = Jijk, (85) 



which defines the connection V in terms of 7 and g by equations 7 1 and 72 The action functional for 7 
and g then reduces to 

F^,u{9,V) = y|^(i?-2A-g|d7|2-||7|2)dy (86) 

= j^iR-'2A-C3\d^\l.form-C4h\l-forJdV (87) 

where we note that the norm of a fully antisymmetric /c-tensor lo is kl times that of the definition of the 
norm of w as a A;-form. Next let 

7 = *{v*) (88) 

where v is a vector field, v* is the 1 form dual to v, and * is the Hodge star operator which takes k forms 
to (4 — A;) forms in dimension 4. We note that |7l3_jorm ~ ^'^'^ Id-lll-form ~ ' '")'^ because 
of the Lorentzian metric. With the above substitution, the action functional becomes 

F^,u{9, V) = / {R-2K + C3(V • vf + cM^) dV, (89) 
Ju 
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where V • v denotes the divergence of v. This amounts to a change of variables for the action from {g, 7) 

to {g,v). 

Computing the Euler-Lagrange equation for variations of the vector field v is straightforward using 
the divergence theorem, but computing the Euler-Lagrange equation for variations in the metric g re- 

^( 



quires the following formulas. Let g{s) be a family of metrics, and let 4-g = h. Then direct calculation 



shows that 

= -{Ric,h) + V -{V ■h-V{{g,h))) (90) 

ds 

^dV = l{g,h)dV (91) 

ds 2 

d 1 

-V-v = -{v,V{g,h)) (92) 

ds 2 

-f- |vp = {v®v,h) (93) 
ds 

where v = v* is the 1-form dual to the vector field v. It follows that the Euler-Lagrange equations for 
the action in equation [89] are 



G + Ag = CiV®v -]^{c^{V -vf + C4\v\^) g (94) 

V{V.v) = —V. (95) 

C3 

For the dominant energy condition to be satisfied, we need 03,04 > 0. To arrive at a nontrivial 
equation for v we need C3 7^ and to arrive at a deterministic equation for v we need C4 7^ 0. Hence, 
let's take 03,04 > 0. In this case, these equations are more easily understood if we define a new function 
/ such that 

/ \ 1/2 



so that by equation 95 



\ 1/2 

- V/. (97) 

C4/ 

Using the above two equations and substituting into the original Euler-Lagrange system of equations, we 
get that 

G + Kg = C3 + (98) 

□/ = ^/ (99) 

C3 

has a solution if and only if the original system does. Note that □ is the Laplacian with respect to the 
Lorentzian metric g so that, modulo reparameterizations of the metr ic, the above equations are hyperbolic 
for the scalar function / and the metric g. We note that the above system has its own action 

h,u{9, V)= [ {R-2A- C3\dff - c^f) dV. (100) 
Ju 

Note that F is not quite the same action as F since the last two terms have opposite signs. 
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By equations 7 1 72 85 88 and 97 the connection will have components 



ijk 



1/2 



{*df)ijk + - (gikj + gjk,i - 9ij,k) 



(101) 



in this representative case. We comment that if we had chosen a different quadratic expression of D 
before "with cross terms" between the fully antisymmetric 7 part of D and the trace parts of D, this 
formula would be modified. However, in all of these cases, we still get the Einstein-Klein-Gordon 
equations with a cosmological constant being the relevant Euler-Lagrange equations. The only difference 



is the interpretation of what the connection is in equation 101 However, so far it is not clear how the 
connection manifests itself physically, other than gravitationally, which is effectively the same in all of 
these cases. 

As a last step, we introduce new constants T and fiQ defined such that 



C4 
C3 



and C4 = IQirfiQ 



(102) 



Plugging these new constants into equations 98 and 99 give us the Einstein-Klein-Gordon equations with 
a cosmological constant in geometrized units with the gravitational constant and the speed of light set to 
one, which are 



G + Ag 
□/ 



TV 



Y2 



^f\g 



(103) 
(104) 



where G is the Einstein curvature tensor, / is the scalar field representing dark matter, A is the cosmolog- 
ical constant, and T is a new fundamental constant of nature whose value has yet to be determined. The 
other constant /xq is not a fundamental constant of nature as it can easily be absorbed into /, but simply 
is present for convenience and represents the energy density of an oscillating scalar field of magnitude 
one which is solely a function of t. It is perfectly fine to set ^uq equal to one, just as we have done to the 
speed of light and the gravitational constant. Also, note that since 



□/ = V-V/ 



1 



(105) 



equation 104 is hyperbolic in / when the metric has signature ( — h ++). 



C A Bigger Geometric Theory 

There are compelling geometric and physical motivations to generalize Axiom [T] as follows. 

Axiom 2 For all coordinate charts ^ : Q C N ^ and open sets U whose closure is compact and in 
the interior ofQ, {g,V) is a critical point of the functional 



F^,uig,y) 



QuadMuciM' UMUC'UC) dVn4 



(106) 



with respect to smooth variations of the metric and connection compactly supported in U, for some fixed 
quadratic functional QuadMuC with coefficients in M L) C. 
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We remind the reader that 



M = {gij} and C = {Ti,k} and M' = {gi,,k} and C = {Tijk,i} (107) 

as before. 

The geometric motivation for this generalization is that Axiom [T] actually changes completely if 
we redefine C = {^ij^} and C similarly. That is, we get different theories resulting from the axiom 
depending on which of the tensor indices (the first and the third) of F we raise or lower. The theory with 
all lowered indices seems most interesting, so that has been the topic of this paper. The other theories are 
interesting too, but result in a slight modification of the Einstein-Maxwell equations with a cosmological 
constant, which is less exciting than a geometric theory of dark matter. Of course these other theories 
with raised indices are interesting to study as well. 

On the other hand, if the quadratic expression Quadjvf uc allowed to have coefficients in C as well 
as M, then this issue of raising and lowering indices becomes irrelevant as the same theory is produced 
in all cases. Furthermore, this generalized axiom also generalizes these other theories into one bigger 
theory. 

Hence, it is very natural geometrically to prefer this second axiom. We would have introduced this 
axiom earlier, but we did not want to complicate things unnecessarily. Of course everything we have 
already studied about scalar field dark matter will still be a part of this generalized theory, but there will 
be (at least) two additional matter fields as well, which is a physical motivation for being intrigued by 
this latest axiom. 

We observe that actions of the form 

Fcs>,u{g, V) = / {cR - 2A - Quad (VT) - Quad (D)) dV, (108) 
Ju 

up to an irrelevant boundary term, are allowed by these axioms, and we conjecture that these are all of 
the actions compatible with Axiom[2] We leave this as another interesting conjecture to study. 

If the above conjecture is correct, then the resulting matter fields generically will be the irreducible 
components of the torsion tensor T. The fully antisymmetric part, studied in this paper, leads to a scalar 
field model of dark matter. The trace part, which will be a 1 form, is reminiscent of the vector potential 1 
form A of electromagnetism. The final matter field, which is everything left over, can be thought of as a 
3 tensor which is antisymmetric in the first two indices, antisymmetrizes to zero, and has zero trace. We 
have no idea how to interpret this physically, but this is certainly an interesting direction to pursue. 

D The Coldness of Scalar Field Dark Matter in the Homogeneous and 
Isotropic Case 

In this appendix we prove the following theorem, stated originally in section [3] 

Theorem 2 Suppose that the spacetime metric is both homogeneous and isotropic, and hence is the 
Friedmann-Lemaitre-Robertson-Walker metric —dt^ + a{t)'^ds'^, where ds\ is the constant curvature 
metric of curvature k. If f{t,x) is a real-valued solution to the Klein-Gordon equation (equation^with 
mass term T) with a stress-energy tensor which is isotropic, then f is solely a function oft. Furthermore, 
if we let H{t) = a'{t) /a{t) be the Hubble constant (which of course is actually a function oft), and p{t) 
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and P{t) be the energy density and pressure of the scalar field at each point, then 

P € 



where 



p 1 + e 



4T2 



(109) 



(110) 



and 

fb 1 rb 



pit) dt, P = / P(t) dt, H' = , , (111) 



where a, b are two zeros of f (for example, two consecutive zeros). 

Proof: We first prove that / is only a function of t. By equation |3j the stress-energy tensor T of the 
real scalar field is given by 

^/ Jf^df f\df\\A 

T/^o = 2— (^Y2" + / ) 9- (112) 

For this tensor to be isotropic, the first term in the above expression must be isotropic. Hence, df must 
be a multiple of dt. Integrating df along the level sets of t proves that / is constant on the level sets of t. 
Hence, / is a function of t (which proves that the stress-energy tensor of the scalar field is homogeneous 
as well). 

Then since df = f'{t)dt and hence jd/p = -f'{tf, we get 

pW//"o = /W' + ^ and P(t)//,o = + (113) 

where p = T{dt, dt) is the energy density and P = T(e, e) is the pressure of the scalar field, where e is 
any unit vector perpendicular to the time direction. 

By equation |4] and the formula for the metric Laplacian, 

T2/(i) = (114) 
= \9\-^'^d^(\g\'l^ g^^d.f) (115) 

= a{t)-^dt{-a{t)^f{t)) (116) 
= -f"{t)-?,H{t)f{t) (117) 

since 1^1 = |det(ff)| = a{t)^ and H{t) = a'{t)/a{t). We can then use this result in the calculation 
below, 

= t-mnfiit (118) 



T2 / T2 



jj{t) [f{t) + ^H{t)nt)^ dt (119) 
' - 2Y2^'^*0 ^^^^ '^^ ^^^^^ 

'm'dt)(i-^^) (121) 
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where we have integrated by parts twice and have not picked up any boundary terms since /(a 
f{b). It follows that 

u „ / rb 

\'2 



3H' 

and 



■p=( [ f{tfdt](2-^] (122) 

J a 

proving the theorem. 

We comment that the complex scalar field case is very similar to the above statement since a complex 
scalar field with quadratic potential is precisely equivalent to two independent real scalar fields (given by 
the real and complex parts of the complex scalar field). Complex scalar fields can be defined which have 
much less oscillations of the pressure (which occurs when the real and complex parts of the field are half 
a period out of phase with one another with the same magnitudes), but since real fields are special cases 
of complex fields, for example, this phenomenon is not generically the case. 

We also comment about out assumption that the stress-energy tensor of the real scalar field is 
isotropic. If this scalar field were the only matter in the universe, then it would follow from the Einstein 
equation that the stress energy tensor of the scalar field was both homogeneous and isotropic, assuming 
the spacetime was both of these. However, we want our theorem to be compatible with other matter 
as well. In this case, as long as the other matter fields have isotropic stress-energy tensors, it still fol- 
lows from the Einstein equation that the stress-energy tensor of the scalar field dark matter is isotropic. 
Specifically, since dark energy due to the cosmological constant is always isotropic (and homogeneous), 
it follows that a model universe with just dark matter and dark energy which was homogeneous and 
isotropic would imply that the scalar field dark matter's stress-energy tensor was isotropic. Hence, this 
is a very general assumption when modeUng homogeneous, isotropic spacetimes. 
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